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During the last 15 years or so, and since the pioneering work of E. Wigner, 
F. Dyson and M. L. Mehta, random matrix theory, combinatorial and per- 
colation questions have merged into a very lively area of research, producing 
an outburst of ideas, techniques and connections; in particular, this area 
contains a number of strikingly beautiful gems. The purpose of these five 
Montreal lectures is to present some of these gems in an elementary way, to 
develop some of the basic tools and to show the interplay between these top- 
ics. These lectures were written to be elementary, informal and reasonably 
self-contained and are aimed at researchers wishing to learn this vast and 
beautiful subject. My purpose was to explain these topics at an early stage. 
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rather than give the most general formulation. Throughout, my attitude has 
been to give what is strictly necessary to understand the subject. I have tried 
to provide the reader with plenty of references, although I may and probably 
will have forgotten some of them; if so, my apologies! 

As we now know, random matrix theory has reached maturity and oc- 
cupies a prominent place in mathematics, being at the crossroads of many 
subjects: number theory (zeroes of the Riemann zeta functions), integrable 
systems, asymptotics of orthogonal polynomials, infinite-dimensional diffu- 
sions, communication technology, financial mathematics, just to name a few. 
Almost 200 years ago A. Quetelet tried to establish universality of the nor- 
mal distribution (mostly by empirical means). Here we are, trying to prove 
universality of the many beautiful statistical distributions which come up in 
random matrix theory and which slowly will find their way in everyday' s life. 

This set of five lectures were given during the first week of a random ma- 
trix 2005-summer school at the "Centre de Recherches Mathematiques" in 
Montreal; about half of them are devoted to combinatorial models, whereas 
the remaining ones deal with related random matrix subjects. They have 
grown from another set of ten lectures I gave at Leeds (2002 London Math- 
ematical Society Annual Lectures), and semester courses or lecture series at 
Brandeis University, at the University of California (Miller Institute, Berke- 
ley), at the Universiteit Leuven (Prancqui chair, KULeuven) and at the Uni- 
versite de Louvain (UCLouvain). 

I would like to thank many friends, colleagues and graduate students in 
the audience(s), who have contributed to these lectures, through their com- 
ments, remarks, questions, etc... , especially Mark Adler, Ira Gessel, Alberto 
Griinbaum, Luc Haine, Arno Kuijlaars, Walter Van Assche, Pol Vanhaecke, 
and also Jonathan Delepine, Didier Vanderstichelen, Tom Claeys, Maurice 
Duits, Maarten Vanlessen, Aminul Huq, Dong Wang and many others.... 

Last, but not least, 1 would like to thank John Harnad for creating such 
a stimulating (and friendly) environment during the 2005-event on "random 
matrices" at Montreal. Finally, I would label it a success if this set of lectures 
motivated a few young people to enter this exciting subject. 
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1 Permutations, words, generalized permu- 
tations and percolation 

1.1 Longest increasing subsequences in permutations, 
words and generalized permutations 

(i) Permutations tt := 7r„ of 1, . . . , n are given by 



Sn ^ '^n = [ ^ ^ ) ) 1 < Ji, • ■ ■ , in ^ ^ all distinct integers 

\ Jl ■ ■ ■ Jn J 

with TTnik) = jk- Then 

= nl 

An increasing subsequence of 7r„ e S'^ is a sequence 1 < ii < ... < < n, 
such that 7r„(ii) < ... < 7r„(ifc). Define 

L{TTn) = length of a longest (strictly) increasing subsequence of 7r„ . 

(1.1.1) 

Notice that there may be many longest (strictly) increasing subsequences! 
Question: Given uniform probability on compute 

P-di.^) <k,n,e S„) = #{^('^-.) < ... g S„] ^ , 

(Ulam's problem 1961) 

Example: for Try = ^ ^ ^ ^ ^ g y 5 ^ ' have ^(717) = 4. A longest 
increasing sequence is underlined; it is not necessary unique. 



(ii) Words tt := tt^ of length n from an alphabet 1, . . . ,q are given by 
integers 

\ Jl ■ ■ ■ Jn J 



with 7r«(A;) = jk- Then 
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An increasing subsequence of vr.^^ G S"^ is given by a sequence 1 < ii < . . . < 
ik < such that vr^(ii) < ... < vr^X^fc)- before, define 

L(7r^) = length of the longest weakly increasing subsequence of vr^ . 

(1.1.2) 

Question: Given uniform probability on S*^, compute 



Example: for7r=^2 ^ ^ g have L{t{) = 3. A longest 

increasing sequence is underlined. 

(iii) Generalized permutations vr := tt^''' are defined by an array of 
integers 

ii ... ir 



n ^ n 



Jl • • • 3n 



subjected to 

l<ii<i2< ■■■<in<P and 1 < ji, . . . , jn < g 
with 2fc = ik+i implying jk < jk+i- 

Then 

#GPr = ( - ' 

An increasing subsequence of a generalized permutation vr is defined as 



, . . . , 



C TT 



with Ti < . . . < rm and jV-i < Jr2 ^ • • • ^ jr^ ■ Define 

L{tt) := length of the longest weakly increasing subsequence of vr. 



Example: For(;^;^;^^^^':|'^'^^)G GF^q, we have L{it) = 5. 



1112223344 
2331221213 

For more information on these matters, see Stanley [80 
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1.2 Young diagrams and Schur polynomials 



Standard references to this subject are MacDonald[68], Sagan[77j, Stanley [80]. 
Stanton and White [81]. To set the notation, we remind the reader of a few 
basic facts. 



A partition X of n (noted A h n) or a Young diagram A of weight n is 
represented by a sequence of integers Ai>A2>...>A£>0, such that 
n = |A| := Ai + ... + Xi, n = \X\ is called the weight. A dual Young 
diagram X^ = ((A^)i > (A^)2 > ...) is the diagram obtained by flipping 
the diagram A about its diagonal; set 



A 



length of ith column of A. 



:i.2.i) 



Clearly |A| = |A^|. For future use, introduce the following notation: 



Y 

YP 
YS 



{all partitions A} 

{all partitions X\- n} 

{all partitions, with A^ < p} 

{all partitions A h n, with A^ < p} 



[1.2.2) 



• A semi-standard Young tableau of shape A is an array of integers > 
placed in box {i,j) in the Young diagram A, which are weakly increasing 
from left to right and strictly increasing from top to bottom. 

• A standard Young tableau of shape A h n is an array of integers 1, ...,n = 
\X\ placed in the Young diagram, which are strictly increasing from left 
to right and from top to bottom. For A h n, define 

n\ _ n { standard tableaux of shape A 1 
■~ ^ I filled with integers 1, |A| /' 

• The Schur function sa associated with a Young diagram A h n is a 
symmetric function in the variables Xi,X2i (finite or infinite), defined 
as 

~ / N -TT #{times i appears in P} /. r, o\ 

SA(xi,a;2,...) = 2^ VV^i (^-^-S) 

semi- standard % 
tableaux P of shape A 
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It equals a polynomial Sx{ti, t2, ■ ■ •) (which will be denoted without the 
tilde) in the symmetric variables ktk = X]j>i^i^5 

sx{xi,X2, ...) = sx{ti, t2, . . .) = det (sA.-»+j(t))i<ij<^ , (1.2.4) 

for any m > n. In this formula Si{t) = for i < and Si{t) for z > is 
defined as 

eEr^^ :=J2s.^h,t2,...)z\ 

Note for X\- n, 

Sx{Xi, X2, ...)= Xi ... X„ + .... = 1^ + • • • 

Given two partitions A ^ /i, (i.e. Aj > fii), the diagram A\/i denotes the 
diagram obtained by removing form A. The skew-Schur polynomial 
Sx\^ associated with a Young diagram A\/i h n is a symmetric function 
in the variables Xi,X2, (finite or infinite), defined by 

#{times i appears in P} 

X; 



SAv(a;i,X2,...) = XI n 

semi- standard i 
skew- tableaux P 
of shape X\fi 

= det(sA,_^^._.+,(t))^<^_^.<„. 
Similarly, for A\/i h n, 

Sx\i,{Xi,X2, ■■■) = /^^^ Xi . . .Xn + . . . . 

The hook length of the i, jth box is defined h^j := Aj + Aj — i — j + 1. 
Also define 



n 

A.m(m + Ai - 1, . . . , m + Am - m) 



rm/ 

rn -i- A ■ — 7 11 

for m > A7.(1.2.5) 



^h^j := Xi + Xj — i — j + 1 is the hook length of the i, jth box m the Young diagram; i.e., 
the number of boxes covered by the hook formed by drawing a horizontal hne emanating 
from the center of the box to the right and a vertical line emanating from the center of 
the box to the bottom of the diagram. 
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The number of standard Young tableaux of a given shape A = (Ai > 
... > A„) is given b}0 



standard tableaux of shape A 
filled with the integers 1, . . . , |A| 

coefficient of xi . . . x„ in sx(x) 



|A|! . . . 



lAlldet 



^sa(u,0,0, ...) 



Ai-1,. 

m 

Y[{m + A, 



l<i j'<?n 



Am(m + Ai - 1, . . . , m + Am - m) -p 
|A|! — , for any m > A^ . 



- I 



(1.2.6) 

In particular, for any m > Xj and arbitrary m G M, 



s,{u, 0, 0, . . .) = = JA| Am(m + A,-l,...,m + Am-m 



|A|! + 

(1.2.7) 



The number of semi-standard Young tableaux of a given shape A h n, 
filled with numbers 1 to g for g > 1: 



# 



semi-standard tableaux of shape A 
filled with numbers from 1 to g 
<? 

Sa(C^,0,0,...) 



n 



j - i + g 



^using the Vandermonde determinant Am{zi, . . . , z„i) ~ ni<i<j<m(^2 ~ ^j) 
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( Ag{q + \i-l,...,q + X,-q) 



n 

i=l 



when q > Xj , 



0, when q < Xj 



using the fact that 



n (-^^ - ^ + ^) 



nLi(g + A.-^)! 



:i.2.8) 



:i.2.9) 



Pieri's formula: given an integer r > and the Schur polynomial s^, 
the following holds 



SASr 



[1.2.10) 



/^\A— horizontal strip 
and |/^\A|— r 



Note /i\A is an horizontal r-strip, when they interlace, /ii > Ai > > 
. . ., and |/i\A| = r. 



1.3 Robinson- Schensted-Knuth correspondence for gen- 
eralized permutations 

Define the set of p x g integer matrices (see [SD] ) 



Mat^'"^ := Iw = {Wij)i<r<v, E Z>o and 



Wij = n 



Theorem 1.1 There is a 1-1 correspondence between the following three 
sets: 



npp,q 



TT < > 



two semi- standard Young tableaux {P,Q), 
of same, but arbitrary shape X\- n, filled 
resp. with integers (1, . . . ,p) and (1, . . . , g) 



{P,Q) 



W{n) = (wij) i<,<p 



1 < j < g 
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where 



Wij = # times that 



J 



I 



) 



Therefore, we hav^ 
f pq + n — 1 



) 



(1.3.1) 



Also, we have equality between the length of the longest weakly increasing 
subsequence of the generalized permutation it, the length of the first row of 
the associated Young diagram and the weight of the optimal path: 



the correspondence constructs two semi-standard Young tableaux P, Q hav- 
ing the same shape A. This construction is inductive. Namely, having ob- 
tained two equally shaped Young diagrams Pk, Qk from 



with the numbers (ji, ...,jk) in the boxes of Pk and the numbers (ii, ...,ik) in 
the boxes of Qk, one forms a new diagram Pk+i, by creating a new box in the 
first row of P, containing the next number jk+i , according to the following 
rule: 

(i) if jk+i > all numbers appearing in the first row of Pk-, then one creates a 
new box with jk+i in that box to the right of the first row, 




Sketch of Proof: Given a generalized permutation 





p 



^Use the notation sa(1P) — sa(1, • • • , 1, 0, 0, . . .). 
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(ii) if not, place jk+i in the box (of the first row) containing the smallest 
integer > jk+i- The integer, which was in that box, then gets pushed 
down to the second row of Pk according to the rule (i) or (ii), and the 
process starts afresh at the second row. 

The diagram Q is a bookkeeping device; namely, add a box (with the number 
ik+i in it) to Qk exactly at the place, where the new box has been added to 
Pk- This produces a new diagram Qk+i of same shape as Pk+i- The inverse 
of this map is constructed by reversing the steps above. 

Formula ffTXTD follows from ffTXg]) . ■ 

Example: For n = 10, p = 4, g = 3, 



1112223344 
2 3 3 1 2 2 1 2 1 3 



with Lin) 



t 



iP,Q) 



\ 



1 


1 


1 


2|3| 


1 


1 


1 


3|4 


2 


2 


2 




2 


2 


2 




3 


3 






3 


4 







L(7t) 



Ai 



t 



W 



(® - 




2\ 




i 




1 


® 







i 




1 


© - 


- ® 






i 


\ 1 





©/ 



with L(7t) 



E 

(ij)e{path} 



Wij = 5 



The RSK algorithm proceeds as follows: 
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adding \ ^ 2 



2\ f 2 \ f 2\ /3 



2 VI 



133 123 122 

''"^"^ 2 2 3 2 3 3 



1 1 2 
2 
3 



Q 111 m m m 222 

^^^^ 2 22 222 



1 1 1 
2 
3 



(0 (0 (3 

1122 1112 11123 

123 122 122 - { ^ 

3 33 33 ~ \ Q 

1113 1113 11134 

222 222 222 

3 3 4 3 4 

yielding the set (P, Q) of semi-standard Young tableaux above. 
1.4 The Cauchy identity 

Theorem 1.2 Using the customary change of variables ^jt>i -^fc ~ S^>i 
isi, we have 
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Proof: On the one hand, to every vr G GP = Unpg^-^n''^; associate a 
monomial, as follows 



j 



Y\ XiUj (with multiplicities) (1-4-1) 



Therefore, taking into account the multiplicity of ^ J G tt. 



One must think of the product on the right hand side of fll.4.21) in a definite 
order, as follows: 

n^= 

.-^^^ 1 - XiUj 1 - xiyi 1 - xiy2 1 - Si^/s 

1 1 1 



1 - X21/1 1 - X2y2 1 - 3:2^3 



and similarly for the expanded version. Expanding out all the products, 

+ a;i2/, + x\y] + x\y] + . . .) + X22/,- + x^^+x^ + ...)•••, (1-4.3) 

leads to a sum of monomials, each of which can be interpreted as a gener- 
alized permutation, upon respecting the prescribed order. Vice-versa each 
generalized permutation can be found among those monomials. As an exam- 
ple illustrating identity fll.4.21) . the monomial Xiy2x\yl Xg^l X2,yz x^yi X4?/2, 
appearing in the expanded version of fll.4.3p . maps into the generalized per- 
mutation 

'1 1 1 3 3 3 4 4 
2 3 3 2 2 3 1 2 
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and vice-versa. On the other hand, to every couple of semi-standard Young 
tableaux {P,Q), we associate 



^ J~J ^#timeB i appears in Q ^ 



appears in Q | j ^times j appears in P 
j 



Therefore, the Robinson-Schensted-Knuth construction, mapping the gener- 
alized permutation tt into two semi-standard Young tableaux (P, Q) of same 
shape A, implies 



mcs j appears in P 



2 3 



Then, summing up over all vr G GP, using the fact that RSK is a bijection, 
and using the definition of Schur polynomials, one computes 

n ^^^j = 

TreGP 

3 J 

mes j appears in P 



El I ^#tiines i appears in Q I I ^.#ti 
11 lli'j 

all (P,Q)with i j 

shape P— shape Q 

E\ ^ I I #times i appears in Q I I #ti 

11 ^» lli'j 



mes j appears in P 



X all (P,Q)with i j 

shape P=shape Q=\ 



= E 



^ ^ XI * appears in Q ^ ^ J~J ^#ti 



all Q with 



;imes j appears in P 



shape Q = \ / \ shape P = X 



all P with 



^SA(a;)sA(|/), 
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using the definition fll.2.31) of tlie Scliur polynomial. The proof is finished by 
observing that 



.k>l £>1 

t 

I I 
i>l i>l 



I 

k/>l t>l 

= log n ^ ^kVe)'^, 
k,e>i 

ending the proof of Theorem ll.2[ ■ 

1.5 Uniform Probability on Permutations, Plancherel 
Measure and Random Walks 

1.5.1 Plancherel measure 

In this section, one needs 



Mat;^'"(0, l):= {W = (m 



ij)^<i,j<n-i 



with exactly one 1 

in each row and column 

and otherwise all zeros 



See [80], [251 ES] and references within. 

Proposition 1.3 For permutations, we have a 1-1 correspondence between 

two standard Young tableaux {P,Q), 
Sn { of same shape X and size n, ^ ■^=r- Mat"'"'(0, 1). 
each filled with numbers (1, . . . , 



TTn < ^ (-P, Q) < > W{tt) = {Wij)ij>i, 

(1.5.1) 
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Uniform probability on Sn induces a probability P" (Plancherel measure) 
on Young diagrams Y^, given by (m :— Xj) 

( permutations^ i'f^'|2 
P'* (A) = — # <^ leading to } ^ 



shape A J 



n\ \ V \ I nl 

nl Sa(1,0,...)' 



and so 



Am{m + Ai - 1, . . . , m + A^ - m)^ 
= nl 

Yl{m + Xi-i)\ 

(1.5.2) 

Xt-n 

Finally, the length of the longest increasing subsequence in permutation iTn, 
the length of the first row of the partition X and the weight of the optimal 
path L{W) in the percolation matrix W{7r) are all equal: 

{over right/down 
y Wij, paths starting from 
< ^ / \ / \ 

entry (1, 1) to [n, n) 

Hence 

P-{L{n)<e)^ ^ ^-^-n\ Yl SA(l,0,...f. 

A e Y„ ■ A e Y„ 
Xi<e Xi<e 

Proof: A permutation is a generalized permutation, but with integers ii, . . . , i„ 
and ji, ■ ■ ■ ,jn all distinct and thus both tableaux P and Q arc standard. 

Consider now the uniform probability P" on permutations in S^, from the 
RSK correspondence we have the one-to-one correspondence, given a fixed 
partition A, 

standard tableaux 1 ( standard tableaux 



permutations 
leading to 
the shape A 



leading to > ^ < ^^^P^ ^' [xl ^^^P*" ^ 

' ' filled with integers | | filled with integers 

l,...,n J [ l,...,n 
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and thus, using formulae fll.2.6p and fll.2.8p and noticing that 8^(1"^) = for 
A7 > q, 



pn (\\ _ ^ _LL j permutations leading 1 _ {f^Y 
^ ' " 1 to the shape A J 



AG 



with 



agy„ 



Formula fll.5.2p follows immediately from the explicit values (I1.2.6P and 
ffT:2:8D for ■ 



Example: For permutation tts 
rithm gives 

P =^ 5 1 



1 2 3 4 5 
5 14 3 2 



G ^5, the RSK algo- 



Hence 



iP,Q) 





1 


2]\ 


/ 


1 


3 \ 




3 




2 






4 






4 




v 


5 


/ 


V 


5 


/ 



1 


4 


1 3 




1 


2 




5 




4 




3 










5 




4 














5 






1 


3 


1 3 




1 


3 




2 




2 




2 










4 




4 














5 








/© 













1 \ 




i 














® 


















i 














© - 


- ® - 


® 




® - 


- ® 














i 










1 







® 














i 




U 


1 










® / 



Remark: The Robinson-Schensted-Knuth correspondence has the following 
properties 



TT ^ (P, Q), then ^ {Q, P) 
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length (longest increasing subsequence of vr) = # (columns in P) 
length (longest decreasing subsequence of vr) = # (rows in P) 
I, then vr ^ (P, P) 



TC 



I, with k fixed points, then P has exactly k columns of odd 



length. 
1.5.2 Solitaire game 



:i.5.3) 



With Aldous and Diaconis consider a deck of cards 1, . . . ,n thoroughly 
shuffled and put those cards one at a time into piles, as follows: 

(1) a low card may be placed on a higher card, or can be put into a new 
pile to the right of the existing pile. 

(2) only the top card of the pile is seen. If the card which turns up is higher 
than the card showing on the table, then start with that card a new 
pile to the right of the others. 

Question: What is the optimal strategy which minimizes the number of piles? 
Answer , put the next card always on the leftmost possible pile! 

Example: Consider a deck of 7 cards, appearing in the order 3, 1, 4, 2, 6, 7, 
5. The optimal strategy is as follows: 



3 
1 



3 


4 


1 





3 


4 


h- ' 


2 



3 


4 


6 


h- ' 


2 





3 


4 


6 


7 


h- ' 


2 





3 


4 


6 


7 


I— ' 


2 


5 





This optimal strategy leads to 4 piles! For a deck of 52 cards, you will 
find in the average between 10-13 piles and having 9 piles or less occurs 
approximately 5% of the times. It turns out that, for a given permutation. 



number of piles = length of longest increasing sequence. 
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1.5.3 Anticipating large n asymptotics 

Anticipating the results in section 14.21 and 19.21 given tlie permutations of 
(1, . . . , n), given a percolation matrix of size n and given a card game of size 
n, the numbers fluctuate about 2^/n like 

L{nn) = length of the longest increasing subsequence 

= weight of the optimal path in the (0, l)-percolation matrix 
= number of piles in the solitaire game 
^ 2v^ + n^/^J^, 

where is a probability distribution, the Tracy- Widom distribution, with 

E{J^) = -1.77109 and a(J^) = 0.9018. 
In particular for n = 52 cards, 

E{L{7i52)) ^ 2^52 + (52)i/''(-l. 77109) = 11.0005. 
The Tracy- Widom distribution will be discussed extensively in section 19. 2[ 

1.5.4 A transition probability and Plancherel measure 

Proposition 1.4 ( f9l[l25^ ) Pn on can be constructed from Pn-i on Yn-i, 
by means of a transition probability, as follows 

Pnifi)= ^n-l(AMA,^), fieYn 

AeY„_i 

where 

if X & and /i G Y„ are 

such that /i = A + □ 
otherwise 

^ p(A,Ai) = 1, for fixed X. 



p{X,fi) :-- 





is a transition probability, i.e. 
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Proof: Indeed, for fixed /i, one computes 



AeY„_i AGY„_i V J \J 

such that /J.— A + DgYtt, 

sucl" 



7 f 

such that ^t — A + D^Ytt, 



n! 

Indeed, given a standard tableau of shape A, filled with 1, . . . , ra — 1, adjoining 
a box to A such as to form a Young diagram fi and putting n in that box 
yield a new standard tableau (of shape fi). 

That p{X, fi) is a transition probability follows from Pieri's formula fll.2.10p . 
applied to r = 1, upon putting ti = Sif. 

f'^° 1 

/Je¥„ with XeY„_i with ' ' AGY„_i with 



|A + n|/^ 

^=^+0 A+n6¥„ A+DeY 





|A|! 


\x + n\ 





Corollary 1.5 The following probability 

P„(/il < Xi, . . . < Xfc) 

decreases, when n increases. 
Proof: Indeed 

PnifJ'l <Xi,...,fIk< Xk) 

= E E ^n-l(A)p(A,/i) 



MSYn AGY„_1 
all fJ'{<x^ 



E P„-i(AMA,^) 
AeY„_i,AteY„ 

such that 
all fli^Xi 



V P„-i(A) = P„-i(Ai < xi, . . . ,Xk < Xk) 

ey;_i 

Lch t 

=A- 

proving Corollary II. 51 



< 

A6l 

such that 
all fJ-i'^Xi 
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1.5.5 Random walks 

Consider m random walkers in Z, starting from distinct points x := {xi < 
. . . < Xm)i such that, at each moment, only one walker moves either one step 
to the left or one step to the right. Notice that m walkers in Z, obeying this 
rule, is tantamount to a random walk in Z™, where at each point the only 
moves are 



with all possible moves equally likely. That is to say the walk has at each 
point 2m possibilities and thus at time T there are (2m) different paths. 
Denote by the probability for such a walk, where x refers to the initial 
condition. Requiring these walks not to intersect turns out to he closely related 
to the problem of longest increasing subsequences in random permutations, as 
is shown in the Proposition below. For skew-partitions, see section 11.21 For 
references, see jSOl [TO]. 

Proposition 1.6 

/ that m walkers in Z, \ 



±ei, . . . , ±e, 



P. 



\ 



reach yi < . . . < y„t 
in T steps, without 
ever intersecting 



J 




(1.5.4) 



where jj, 



V are fixed partitions defined by the points Xi and yi 



k-1 - Xk, Vk = k - I - yk 




m 



1 




m 



T 



Tl + Tr, ^{xi - Vi) =Tl- Tr. 
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In particular, close packing of the walkers at times and T implies 

that m walkers in Z reach 
1, ... ,171 in 2n steps, 



Pi. 



without ever intersecting 
1 /2n 



{2m] 



2n 



n 



/ Xl-n 



Ahn 



{2n)\ H^{nn G Sn such that L{Tin) < ^} 



(2m) 



2n 



(1.5.5) 



1 2 3 4 5 6 



1 2 3 4 5 6 

(Six non-intersecting walkers leaving from and returning to 1, . . . , 6. 



Proof: Step 1: Associate to a given walk a sequence of Tl I/'s and Tr i?'s: 
LRRRLRLLR...R, (1.5.6) 
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thus recording the nature of the move, left or right, at the first instant, at 
the second instant, etc... 

If the k^^ walker is to go from Xk to yk, then 

J right moves 1 J left moves 

yk-Xk = #<^ for k'^ walker J ^ \ for k"" walker 

and so, if 

rp ._ ni Isft moves 1 ^ ._ » f right moves 

^ ' y for all m walkers J ^ I ^ walkers 

we have, since at each instant exactly one walker moves, 

m 

Tr + Tl = T and Tr-Tl = ^(z/fc - a^fc), 

1 

from which 



Next, we show there is a canonical way to map a walk, corresponding to 
(11. 5. 61) into one with left moves only at instants 1,...,Tl and then right 
moves at instants + 1, . . . , + Tr = T, thus corresponding to a sequence 



Tl 



L L L ... L R R R ... R . (1.5.7) 



This map, originally found by Forrester [12] takes on here a slightly different 
(but canonical) form. Indeed, in a typical sequence, as ( ]1.5.6p . 

L R R R L R L L R . . . R , (1.5.8) 

consider the first sequence R L (underlined) you encounter, in reading from 
left to right. It corresponds to one of the following three configurations (in 
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the left column), 



L 
R 



R 
L 



L 
R 



R 
L 



L 
R 



R 
L 



which then can be transformed into a new configuration L R, with same 
beginning and end, thus yielding a new sequence; in the third case the re- 
flection occurs the first place it can. So, by the moves above, the original 
configuration fll.5.6p can be transformed in a new one. In the new sequence, 
pick again the first sequence RL, reading from left to right, and use again 
one of the moves. So, this leads again to a new sequence, etc... 



L 


R 


R 


R 


L 


R 


L 


L 


R . 


. R 


L 


R 


R 


L 


R 


R 


L 


L 


R . 


. R 


L 


R 


L 


R 


R 


R 


L 


L 


R . 


. R 


L 


L 


R 


R 


R 


R 


L 


L 


R . 


. R 


L 


L 


R 


R 


R 


L 


L 


L 


R . 


. R 



L L L ... L R R R ... R ; (1.5.9) 

Since this procedure is invertible, it gives a one-to-one map between all the 
left-right walks corresponding to a given sequence, with Us and Tj^ i?'s 

LRRRLRLLR...R; (1.5.10) 

and all the walks corresponding to 

L L L ... L R R R ... R . (1.5.11) 
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Thus, a walk corresponding to (11. 5. Ill) will map into 




different walks, 



corresponding to the ^^^^^^ j number of permutations of L's and Tr i?'s. 

Step 2: To two standard tableaux P, Q of shape A = (Ai > . . . > Am > 
0) we associate a random walk, going with (11. 5. lip , in the following way. 
Consider the situation where m walkers start at 0,1,. ..,m — 1. 



instants of left move 



1*^* walker 


Cll 


C12 


Cl3 


Ci4 


Cl5 


C16 


2°^^ walker 


C21 


C22 


C23 








3'''^ walker 


C31 


C32 


C33 










C41 













• The walker starts at and moves to the left, only at instants 

Cii = content of box(l,i) G P 
and thus has made, in the end, Ai steps to the left. 

• The k^^ walker {I < k < m) starts at A; — 1 and moves to the left, only at 
instants 

Cfcj = content of box(A;, i) G P. 
and thus has made, in the end, A2 steps to the left. 



• Finally, walker m = walks according to the contents of the last row. 

Since the tableau is standard, filled with the numbers 1, . . . , n the walkers 
never meet and at each moment exactly one walker moves, until instant 
n = |A|, during which they have moved from position 

0<1< ... < k - I < ... < m-l 
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to position 



-Ai + < -A2 + !<...< -Xk + k - 1 < . . . < -Xm + m - 1 

That is to say the final position is given by unfolding the right hand boundary 
of A, the horizontal (fat) segments refer to gaps between the final positions 
and vertical segments refer to contiguous final positions. 

In the same fashion, one associates a similar walk to the other tableau Q, 
with the walkers also moving left. These walkers will have reached the same 
position as in the first case, since the final position only depends on the shape 
of P or Q. Therefore, reversing the time for the second set of walks, one puts 
the two walks together, thus yielding m non-intersecting walkers moving the 
first half of the time to the left and then the second half of the time to the 
right, as in the example below. Therefore the number of ways that m walkers 
start from and return to 0, . . . , m — 1, without ever meeting each other, by 
first moving to the left and then to the right, is exactly \\-n (f^) . 



{P, Q) = 




More generally, an analogous argument shows that the number of ways 
that walkers leave from xi < . . . < Xm and end up at 7/1 < . . . < t/^ at time 
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T, without intersection and by first moving to the left and then to the right, 
is given by 

J2 /'^'^/'^'^ (1-5.12) 

Xj<m 

On the other hand, there are (^^^^^^^ sequences of L's and Tr i?'s, which 

combined with fll.5.12p yields formula (11. 5. 41) . 

In the close packing situation, one has yUfc = = for all k, and so 
/i = z/ = and = = T/2. With these data, (11. 5. 51) is an immediate 
consequence of (11.5.41) . ■ 



1.6 Probability measure on words 

Remember from section 1.1 words vr := vr^ of length n from an alphabet 
1, . . . , g and from section II. 2[ the set of partitions A h n, with < g}. 
Also define the set of n x g matrices. 



Mat;''(0,l):= {W={w,,),<,,< 



with exactly one 1 

in each row 

and otherwise all zeros 



For references, see 

Proposition 1.7 In particular, forwards, we have the 1-1 correspondence 

semi-standard and standard Young tableaux 
{ (P^Q) of same shape and of size n, filled \ •i=?- Mat^'^(0, 1). 
resp., with integers (1, . . . , g) and {1, . . . ,n) 



TT < ^ (P,g) < > W{7r) = {Wij)ij>i. 

(1.6.1) 

Uniform probability P"''' on S'^ induces a probability P"''' on Young diagrams 
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A G Yn \ given by 

( words m SI i /a- /igN 



^ ' to shape A I 

n! , ^ , q q 

= -sa(1,0,...) SA(g,|,|,...) 

n! Ag(g+Ai-l,...,g+Ag-g)^ 

g^n^! n^^+^'~')' 
1 1 

(1.6.2) 

i^Sl = J2f' Sx{l') = q''. (1.6.3) 

Ahn 

Finally, given the correspondence 1^1.6. the length L^ir) of the longest 
weakly increasing subsequence of the word vr equals 

fover right/down 
y Wij, paths starting from 



patfts I en^ry (1,1) to (n, g) 



and i/iits 



P'^'^(L(vr) < £) = 5^ ^s,(l,0,...)s,(g,|,| 

Ai < £ 
A G Yi'^) 



5 • • V 



Proof: A word is a special instance of generalized permutation, where the 
numbers ii, . . . ,in are all distinct. Therefore the RSK correspondence holds 
as before, except that Q becomes a standard tableau; thus, a word maps 
to a pair of arbitrary Young tableaux {P,Q), with P semi-standard and Q 
standard and converse. Also the correspondence with integer matrices is the 
same, with the extra requirement that the matrix contains and I's, with 
each row containing exactly one 1. 
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Consider now the uniform probability P^''^ on words in S'^; from the 
RSK correspondence we have the one-to-one correspondence, given a fixed 
partition A, 



semi- standard 
tableaux of shape 
A, filled with 
integers 1, . . . ,q 

and thus, using formulae (ll.2.6p and (I1.2.8P and noticing that 8^(1'^) = for 
Xj > q, 



words in S'^ 
leading 
to shape A 



standard tableaux 
of shape A 
filled with integers 
1 n 



words leading 
to the shape A 



with 



J2 ^"'"(A) = 1. 



Formula (11. 6. 21) follows immediately from an explicit evaluation f ll.2.6p and 
([L21D for and Sa(1''). ■ 



Example: For word 



TT 



the RSK algorithm gives 

2 1 
2 



1 2 3 4 5 

2 113 2 



1 1 
2 



G S; 



1 3 
2 



Hence 



TC 



1 


1 


2 




1 


3 


4 


2 


3 


1 


2 


5 





1 3 




1 


1 2 






2 


3 


3 4 




1 


3 4 






2 


5 


/© 


1 




o\ 


i 








® 










i 








® - 






- ® 


















® 








i 




1 




® / 
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and L(7r) = Ai = L{W) = 3. 



1.7 Generalized Permutations, Percolation and growth 
models 

The purpose of this section is to show that uniform probabihty, generahzed 
permutations, percolation, queuing and growth models are intimately related. 
Their probabilities are ultimately given by the same formulae. 



1.7.1 Probability on Young diagrams induced by uniform proba- 
bility on generalized permutations 

Proposition 1.8 (Johansson 156]) Uniform probability P^''^ on GP^''^ in- 
duces a probability P^'*^ on Young diagrams Y™"''^'^\ given by 

~„„, , 1 J generalized permutations 



jj^GP^'^ \ leading to shape A 
§a(F)§a(1'^) 

i^GP'n' 

\{q + Ai - 1, g + A, - n M q < V. 



^GP'n'= E §a(F)s,(P)= ('^^+^^ M. (1.7.1) 

Proof: By the RSK correspondence, we have the one-to-one correspondence, 

{semi- standard "1 ( semi-standard 

tableaux of shape I I tableaux of shape 

A, filled with I I ^^^^^ with 

integers 1, . . . , g I I integers 1, . . . , p 



wit, 



* There is no loss of generality in assuming q < p. 



31 



and thus, for A G Y™™^'''^-*, we havf^ 



pp,g^^ = 1 ^ / ^ e GP^ leading 1 ^ sJl'^ls.fF) 

l^J #GPP'^^\ totheshape A j #GP^'^^^^ 

when g < A^, we have automatically sa(1'^) = 0, and thus 



^gYmi„(5,p) 



Notice thal^for m > Xj , 

~ nrr^^ _ + Aj - 1, ...,m + A^ - m) _ -p-j- Aj - Aj + j - 2 



n 



? — ? 

^1 l<i<j<m 



1 



(1.7.2) 



Without loss of generality, assume p > q; then Xi = ii q < i < p. Setting 
hj = Xj + q — j for j = 1, q, we have /ii > ... > /ig = Ag > 0. We now 
compute, using (11.2.81) . 



Sa(F)sa(1 



Xi - Xj + J - l\ YJ TT -^i + J - ^ 



n e^^^i^ nn 

n (/ij - hjf -pi- hj+j-q 
(7-2)2 11 11 



n n - ') n - 

i=l j=q+l l<i<i<<? 

j=0 'i < J J l<i<j<q i=l * 



'^Remember the notation from section 1.2: A G Y"""*-'''^-' means the partition A h n 
satisfies Xj < p, q. 

'Noteni<,<,<,(.7-*)=n|Ioj!- 
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using 



n o-^n n (^-^) 

l<i<j<q i=l i=g+l 

q-l q 

= llJ^■ll{q+l-^){q + 2-^){q + 3- i)...{p - i) 

i=l 

= M + iy-{q + 2)\ {p-iy. 

^' 1! 2! "\p-q-iy '^^ 



9-1 
3=0 

This ends the proof of Proposition 11.81 



1.7.2 Percolation model with geometrically distributed entries 

Consider the ensemble 

Mat^P'^) = {px q matrices M with entries Mjj = 0, 1, 2, . . .} 
with independent and geometrically distributed entries, for fixed < < 1, 
P(M,, = fc) = (1 - Oe', fc = 0,l,2,.... 
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Theorem 1.9 (Johansson l5E]) Then 

{over right/ down 
y Mij, paths starting from 
entry (1, 1) to {p,q) 

has the following distribution, assuming q < p, 

p{L{M)<i) = J2 (i-ere'"§A(i')§A(F) 



hem i=i 

max(h.^ ) <£+g — 1 



where 

q-l 

Zp,, = ^"^(1 - O-'^U^Kp - Q + jY (1-7.3) 

3=0 

Proof: Then the probabihty that M be a given matrix A = (aij) l<^<p equals 

l<j<q 

P i M = {uij) i<i<p j = |T P{Mij = aij), using independence, 

l<j<9 

= n (1 - ^)^"" 



l<i<p 
l<j<q 



Yl l<i<p '^ij 



This probabihty only depends on the total weight \A\ = V. . Oj,-. Hence the 
matriceCl M G MatJ^'' have all equal probability and, in particular, due to 
the fact that, according to Theorem II. ![ the matrices in Mat^*^ are in one-to- 
one correspondence with generalized permutations of size n, with alphabets 
1, . . . ,p and 1, . . . , g, one has 

all an \ 1<J<<3 / 



"^M e Matfj' C Mat^', means that J2 = n. 
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We now compute 



Hence, 

P{L{M) < 



P{L{M) < i \M\ = n) 

#{M G MatP'^ L{M) < 

#{7r e GPr, ^(vr) < i} 

J2 Sa(1^)Sa(F). 



|A|=n 



CXD 



< 



n=0 

oo 



|M| =n) P{\M\ =n) 



E E ^^§A(i^)§A(F)(#GPr)(i-err 



n=0 ^i<f 

I A I —n 



5^ j]s,(i'')s,(F)(i-erei'| 

n=0 Ai<£ 

I A I — n 

5^s,(i'^)s,(F)(i-erei'|. 



Ai<« 



Now, using the expression for sa(1^)sa(1^) in Proposition II. 8[ one computes, 
upon setting hi = q + Xi — i, and noticing that ^ > Ai > . . . Ag > imphes 
i + q-l>h>...>hg>0, 



P{L{M) < 



= J](l-enl'lsA(F)§A(l^) 
AeY9 

Ai<« 
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(1-0 ^^[[jKp-q + jV- 

j=0 

hem i=i 

since the expression in the sum is symmetric in hi, hq. The normahzation 
Zp^q is as announced, ending the proof of Theorem ll.Qi ■ 

1.7.3 Percolation model with exponentially distributed entries 
Theorem 1.10 15^ Consider the ensemble 

Maf''^ = {p X q matrices M with entries} 
with independent and exponentially distributed entries, 

P{Mij <t) = l- e-\ t > and G M. 

Then 

{over right/lower 
y Mij, paths starting from 
' ^ / \ / \ 

entry (1,1) to {p,q) 

has the following distribution, (assuming q < p, without loss of generality), 

/(o,t)9 ^9(-^i' • • • ' -^q) Yli=l -^f '^^ ^dxi 



/(o,oo)9 ^^(^i' • • • ' •^9)^ 11^'=! -^f ''^ ^^dxi 
^ ^ {detMf-'^e-^'^dM, 



P{L{M) < t) 



JMeHq, Sped rum {M)<t 

Remark: It is remarkable that this percolation problem coincides with the 
probability that the spectrum of an Hermitian matrix does not exceed t, 
where the matrix is taken from a (positive definite) random Hermitian en- 
semble with the Laguerre distribution, as appears in the second formula; this 
ensemble will be discussed much later in section I8.2.2[ 
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Proof: For fixed < ^ < 1, let have a geometric distribution 

P{X^ = A;) = (1 - < e < 1, = 0, 1, 2, . . . ; 

then in distribution 

fore^l 

where y is an exponential distributed random variable. Indeed, setting e :- 



0<k<t/s 

£(1 — e)^'^'^ (Riemann sum) 



0<fc<t/£ 



f ds = P{Y < t) 
Jo 



Then, setting ,^ — 1 — £, t = ehi = xt in formula fll.7.3p of Theorem 11.91 
and letting e — > , one computes 

lim P^(L(M) < t) 

e— >0 

= \un 

q /p-g \ 

^ A,(/ii£, . . . , n n + ^) (1 - ^)''''^' 

max(/iie)<fe+((j-l)£ i=l \fe=l / 



^ J\o.t]i Z_T 



(det M)P-'?e-^'"*^dM. 



1 

cctrum (M)<t 

This last identity will be shown later in Proposition I7.2[ ending the proof of 
Theorem 11.101 ■ 
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1.7.4 Queuing problem 

Consider servers 1,2, ... ,q waiting on customers 1, . . . ,p, with a ^^first-in 
first-ouf service rule; see e.g., |2l] and references within. At first, the system 
is empty and then p customers are placed in the first queue. When a customer 
k is done with the i^^ server, he moves to the queue waiting to be served by 
the i + 1*^* server. Let 

V{k,i) = service time of customer k by server i 

be all geometrically independently distributed 

P{V{k,i) = t) = {l~Oe, t = 0,l,2,... . 



servers: 



customers: 



' 1 ' 



Servers 1,2, ... ,q waiting on customers 1, . . . ,p 



Theorem 1.11 The distribution of 



D{p,q) 



departure time for the last 
customer !, at the last server 
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is given by ( assuming q < p) 

PiDip,q)<i) = Yl (l-ere'"s,(l^)s,(P 



heNI i=i ' 



:i.7.4) 



Proof: We show the problem is equivalent to the percolation problem dis- 
cussed in Theorem 11.91 Indeed: 

Step 1: Setting D{p, 0) = -D(0, q) = for all p, q, we have for all p,q > 1, 

D{p, q) = max{D{p - 1, g), D{p, q - 1)) + Vip, q). (1.7.5) 

Indeed, if D{p — l,q) < D{p, q — 1), then customer p — 1 has left server q by 
the time customer p reaches q, so that customer p will not have to queue up 
and will be served immediately. Therefore 

Dip, q) = D{p, g - 1) + V{p, q). 

Now assume 

D{p-l,q)>D{p,q-iy, 

then when customer p reaches server q, then customer p — 1 is still being 
served at queue q. Therefore the departure time of customer p at queue q 

D{p, q) = D{p - 1, g) + V{p, q). 

In particular 

D{p, l)=D(p-l,l) + V(p,l) 

and 

D(1,2) = D(1,1) + K(1,2), 

establishing (11. 7. 51) . 

Step 2: we now prove 



D{p,q) = inax <^ V\/(z,j; 

alt such I t * 
paths ^ 



over right/down paths 
from entry (1,1) to {p,q) 
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where the paths are taken in the random matrix 

V = {V{z,j))^<.<,. 

l<J<q 

By a straightforward calculation, 

D{p,q) = max{D{p-l,q) + V{p,q),D{p,q-l) + V{p,q)) 

= max ( max Vjhj) + V{p, q), 

\(oM'At'^i,<i) path 



,11 paths f ^ 



all paths 
(0,0)-(p,5-l) path 



max 

all paths 
(0,0)^(p,?) 



^path 



ending the proof of Theorem 11.111 ■ 

1.7.5 Discrete polynuclear growth models 

Consider geometric i.i.d. random variables ^^(x, t), with x G Z, t G Z+, 

P(,;(x,t) = fc) = (l-Oe', A;GZ+, 

except 

t) = if t — a; is even or |x| > t. 

Define inductively a growth process, with a height curve h{x,t), with x G 
Z, t G Z+, given by 

/i(x,0) = 0, 



h{x, t + 1) = maxyh{x — 1, t), h{x, t), /i(a; + 1, t)j + a;(x, t + 1). 
For this model, see [59], [79] and [66] . 

Theorem 1.12 The height curve at even sites 2x at times 2t — l is given by 

{over right/down 
y^V{i,j), paths starting from 
entry (1, 1) to [t + x,t — x) 

40 



where 

Thus h{2x, 2t - 



V{hj) ■= ^ii- 3,i+3 - !)• 
1) has again the same distribution as in l[1.7.!^ . 



Proof: It is based on the fact that, setting 

:= h{q-p,q+p- 1), 

one computes 

G{q,p) = max{G{q-l,p),G{q,p-l)) + V{q,p) 

{over all right/down 
Vij, paths starting from 
entry (1, 1) to {q,p) 

So 

h{2x, 2t-l) = G{t + X, t - x), 
establishing Theorem 11.121 



The figure below gives an example of such a growth process and also 
it shows that h{2x,2t + 1) is given by the maximum of right/down paths 
starting at the upper- left corner and going to site {t + x,t — x), where x is 
the running variable along the anti-diagonal. 
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2 Probability on partitions, Toeplitz and Fred- 
holm determinants 



Consider variables x = (xi, X2, ■ ■ ■) and y = {yi, y2, ■ ■ ■), and the correspond- 
ing symmetric functions 

t= {ti,t2,...) and s= (si,S2,...) 

with 



ktk = ^ and ksk = ^ y- 



.k 

I 

i>l i>l 



Following Borodin-Okounkov-Olshanski (see [221 [211 125]), given arbitrary 
X, y, define the (not necessarily positive) probability measure on A G Y, 

P.,yW := ^Ux)sx{y) = ^^x{t)sx{s) (2.0.1) 

with 

Indeed, by Cauchy's identity, 

^Px,yiX) = ^sxix)sxiy)Y[il-Xiyj) 

AeY A l<i,j 

= Y[{l-x,y,)-'Y[{l-x,y,) = l. 

The main objective of this section is to compute 

P.,yiXi<n), (2.0.2) 

which then will be specialized in the next section to specific x's and y's or, 
what is the same, to specific t's and s's. This probability (12.0. 2p has three 
different expressions, one in terms of a determinant of a Toeplitz determinant, 
another in terms of an integral over the unitary group and still another in 
terms of a Fredholm determinant. The Toeplitz representation enables one 
to compute in an effective way the probability (12.0.21) . whereas the Fredholm 
representation is useful, when taking limits for large permutations. In the 
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statement below, we need the Fredholm determinant of a kernel K{i,j), with 
2, j G Z and restricted to [n, oo); it is defined as 

det (l - K{i,j) ] := f^i-ir 52 det{K{x„xj))i<i,,<m. 

\ \n,n+l,...] I 

N ^ 771=0 n<x^<...<xm 

(2.0.3) 

Now, one has the following statement [i5l [23l[26l[85l[8] : 

Theorem 2.1 Given the "probability measure" ^2. 0. the following proba- 
bility has three different expressions: 



\Jsi 2niz 

JU{n) 



l<k,e<n 



det I-K{k,. 



(2.0.4) 



[n,n+l,...]. 



where K{k, i) is a kernel 



Kik, 



1 

27ri 



dz dw e^^^)-^^"') 

|^|=p<l J|^|=p-l>l Z-W 



for k,ieZ 



k-. 



1 

2^1 



l=p<i z^~^^w~^ 

|z|=p-l>l 



dz dw zfVjz) - w£V{w) ^v'(.)-yH 
z — w 

for k,i & Z, with k i, 

(2.0.5) 



with 



Viz) := Vt,siz) := - ^^(t^^"'' - s,z'). 

i>l 



2.1 Probability on partitions expressed as Toeplitz de- 
terminants 

In this subsection, the first part of Theorem 2.1 will be reformulated as 
Proposition 12.21 and also demonstrated: (see e.g. Gessel pS], Tracy- Widom 
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Adler-van Moerbeke [S]) 
Proposition 2.2 Given the "probability measure" 

P(A) = Z-hx{t)sx{s), Z = e^»>i**»^% 
the following holds 

and 

P{X wtth Xl<n)= Z-Uei ( i A^.^-^e^TiUz^+s^^-^) 



l<kl<n 



I SI 27riz J i<fc_^<„ 

Proof: Consider the semi-infinite Toeplitz matrix 

moo(t,s) = (/ifc^)fc,£>o, with s) = / ^'^-^e^rCt^^^-^j^"^) 



27iiz 



Note that 



dtm Js^ 2mz 



OSm isi 2mz 

with initial condition fikei^, 0) = Sm- In matrix notation, this amounts to 
the system of differential equationq^ 

°° = A*moo and °° = — moo(A^)*, with initial condition moo (0,0) = h 

The solution to this initial value problem is given by 
(i) 

moo{t, s) = {fike{t, s))k,e>o (2-1-1) 



^The operator A is the semi-infinite shift matrix, with zeroes everywhere, except for 
I's just above the diagonal, i.e., (Aw)„ = Wn+i- loa is the semi-infinite identity matrix. 
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and 
(ii) 

where 



moo(t, s) = e^"**^'moo(0,0)e 



(2.1.2) 



( 1 Si(t) S2(t) S3(t) ... \ 

1 Si(t) S2(t) ... 

1 si(t) ... 
1 



V ■■ 



(Sj_j^(t^^ l<i<oo . 

1 <j <oo 



Then, by the uniqueness of solutions of ode's, the two solutions coincide, and 
in particular the n x n upper-left blocks of fl2.1.ip and (12.1.21) . namely 



where 



Enit) 



m„(t,s) =E„(t)moo(0,0)E;[(-s) 
/ 1 Si(t) S2(t) ssit) ... s„_i(t) ... \ 

1 Si(t) S2(t) ... S„_2(t) ... 



(2.1.3) 



v 







Si(t) 

1 



(Sj_j(t)) l<i<n . 

l<j<oo 



Therefore the determinants coincide: 

detm„(t,s) = det(E„(t)moo(0,0)E;[(-s)). 



(2.1.4) 



We shall need to expand the right hand side of (12.1.41) in "Fourier series" , 
which is based on the following Lemma: 

Lemma 2.3 Given the semi-infinite initial condition moo(0,0), the expres- 
sion below admits an expansion in Schur polynomials, 

det(E„(t)moo(0,0)ET(-s)) = ^ det(m^'^)sA(t)s,(-s), /orn > 0, 

A, V 

(2.1.5) 

where the sum is taken over all Young diagrams A and v, with first columns 
Xj and vj <n and where 

'^'^ (^A.-j+n,!.,-j+n)i<i j<„. (2.1.6) 



m 
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Proof: The proof of this Lemma is based on the Cauchy-Binet formula, which 
affirms that given two matrices A , B , foT n large > m 

(m,n) (n,m) 

det{AB) = det | ^aubik | 

V « / l<k,e<m 

= det{ak,ii,)i<k,e<mdet{bi^^e)i<k/<m- (2.1-7) 

l<ii<...<im<'n 

Note that every decreasing sequence oo > > ■ ■ ■ > ki > 1 can be mapped 
into a Young diagram Ai > A2 > ... > A„ > 0, by setting kj = j + \n+i-j- 
Relabeling the indices i,j with 1 < i, j < n, by setting j' := n — j + 1, i' : = 
n — 2 + 1, we have 1 < < n and kj—i = Xjr—j' + i' and ki — 1 = Xii—i' + n. 
In other terms, the sequence of integers 00 > A;„ > . . . > /ci > 1 leads to a 
partition Ai = A;„ — n > A2 = kn-i — n + l>...>A„ = A;i — 1>0. The same 
can be done for the sequence 1 < ^1 < ... < ^„ < cxd, leading to the Young 
diagram u, using the same relabeling. Applying the Cauchy-Binet formula 
twice, expression fl2.1.4p leads to: 

det {E^{t) m^{0,0) EJ{-s)) 

J2 det (s.,-.(t)),<,^^.<„det ((moo(0,0)E:(-s)).,,,)^<^_,<„ 

l<fcl<...<fc„<oo 
l<fci<...<fc„<oo 

J2 det(sfc^._i(t))^^^^^.^^ 

l<fei<...<A;„<co 



l<^l<...<^n<00 



J2 det (/ifc^_i,,^,_i)^^^_^.^^det (s,^_,(-s))j 

<...<^n<00 

A6Y 

J2 det (/iA,,_,.+„,,^,,_jv+„)^^^^ ^^^det (s,^,_,,+j.(-s))^^., 



l<i',j'<n 



= ^ det(m^''^)sA(t)s,(-s), 

A, i^gY 

which establishes Lemma [2. 3[ 
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Continuing the proof of Proposition l2.2[ apply now Lemma l^T^ to moo(0, 

detm^'" = det{^x,-i+n,u,-j+n)l<i,j<n ^ 



/oo, leading to 



if and only if 

Xi — i + n = i'i — i + n for all 1 < i < n, 
i.e., A = z/, in which case 

detm^'^ = 1. 
Therefore, from fl2.1.4p . it follows 

y s,it)s,{-s) =det( [ ^z'-'e^^^'^^'-^^^"A 

\Jsi 2mz J ,<fc^^<„ 

So, we have, changing s — s, 

P (A with \l <n) = Z-^ J] sxit)sxis) 

and, using sa(— t) = (— l)l'^ls;^T(t), we also have 
P (A with Ai < n) = ^ sx{t)sx{s) 

= J2 ^X-r{t)SxT{s) 

ASY 



where 

ending the proof of Proposition 12.21 
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2.2 The calculus of infinite wedge spaces 



The material in this section can be found in Kac [61] and Kac-Raina [62] and 
many specific results are due to Borodin-Okounkov-Olshanski (see [23l [2^ 
25]). Given a vector space V = (Bjez^vj with inner-product {vi,Vj) = 6ij, 
the infinite wedge space = A'^V is defined as 

Si > S2 > . . . 

V°° = span ^ f A f S2 A A . . . , 

Sk = —k for A; » 
containing the "vacuum" 

/0 = f_i A t;__2 A . . . . 

The vector space V°° comes equipped with an inner-product ( , ) , making 
the basis vectors Vg^ A f ^2 A . . . orthonormal. To each A; G Z, we associate two 
operations, a wedging with and a contracting operator, removing a Vk, 

^* : ^ ^oo . _ _ _ ^ ^^^^^^^ A . . . A A . . 

i 

Note that 

= 0, if Vk figures in / 

tplivs, A . . .) = 0, if k ^ {si,S2,...). 

Define the shift 

fcez 

acting on V°° as follows 

A^'Vs-, AVs^A... = Vs^+r A A f S3 A . . . 

+Vs-, A Vs2+r AVs-^ A . . . 
+Vs-, A Vs., A Vs,,+r A . . . 

+ ... 
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One checks that 

[A^^fe] = V'.+.,[A^V'fc] = -^^. (2.2.1) 

[A^A^] =£4,-£ (2.2.2) 

and hence 

[J2 UA\ J2 ^^^"1 = -Yl (2.2.3) 

i>l j>l i>l 

Lemma 2.4 (Version of the Cauchy identity) 

gEi>i *i^'e~ '''^ = e^'>i ^ii'^ig- Ej>i '5jA~-''gEi>i uA-^ (2.2.4) 

Proof: When two operators A and B commute with their commutator [A, B], 
then according to Kac [61J, p 308, 

Setting A = ^tiA* and B = -^sjA"^', we find 

i>i j>l 

= g- Ej>i SjA'^ gEi>i *iA'gEi>i 

= pE»>iit.^. p-Y: sjA-^ Y:uA^ 



It is useful to consider the generators of ipi and ipi- 

iJ{z) = J2^'^^, riw) = J2^~'i^j- (2.2.5) 
From fl2XID . it follows that 

k k ^ 

[A^V^*(«;)] = --l^*(«;). 

The two relations above lead to the following, by taking derivatives in tj and 
Si of the left hand side and setting all = Sj = 0: 
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.SI le-^r.- . / -t!J:rtl^!. (2-2.6) 



To each partition A = (Ai > A2 > . . . > A^ > 0), associate a vector 

fx := Vai_i a Va2_2 a ... a WA„-m A A v_(^+2) A . . . e y 

The following lemma holds: 
Lemma 2.5 

eS.>i*^A»/, = ^s^\aW/m, e^»^i*^^-7A = (2-2.7) 



/n particular, 

e^^^^'^^'h = J2^,{t)f, and e^-i^'^^V^ = /g. (2.2.8) 
Proof: First notice that a matrix A e GLoo acts on as follows 



si>si>... 

where 

^si,s2, -- _ f matrix located at the intersection of the 
4,4,... 1^ rows s[, S2, . . . and columns si, S2, . . . of A. 

Here the rows and columns of the bi-infinite matrix are labeled by 
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/ * * * * * \ 



-1 



-2 



-3 V 



Hence, for the bi-infinite matrix e^i° 



S1,S2,. 



det (s,._,J 



Setting Si = Xi — i and defining /Xj in the formula below by s'^ 
one checks 



/a 



(w^, A A . . .) 



E det((< 



s\ >s'>. 



J2 det(s(^^_,)_(A^_,)(t))^^. 

/il-l>/i2-2>... 



j<oo 



The second identity in fl2.2.7p is shown in the same way. Identities (I2.2.8I) 
follow immediately from (12.2.71) , ending the proof of Lemma 12. 5[ ■ 
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We also need: 
Lemma 2.6 For 

= ^'^^ "^a, = ^afciV-i, (2.2.9) 

i>l i>l 

the following identity holds: 

{^a, . . . ^a„*L • • • h) = ((^a,n/0. h)) , , (2-2.10) 

\ / l<k,t<m 

Proof: First one computes for the \l/'s as in fl2.2.9p . 

n^---n.h = E (-i)^^^"^'^det(6.,o,<,,,<^ 

il>...>jm>l 

v^i A . . . Av-i A ... Av-i , A ... A V-i^ A ... A V-i, A ... 

J- ''m I'm — 1 f'Z ''1 

Then acting with the \l/aj. as in fl2.2.9p . it suffices to understand how it acts 
on the wedge products appearing in the expression above, namely: 

'^a. ■ ■ ■'^a v-i A . . . Av-i A ... Av^i , A ... A -O-i, A . . . 

Lti '-I'm -L I'm I'm — 1 ^-l 

= (-l)^^^"^^^det(a.,,),<,,<^^/0 

Thus, combining the two, one finds, using the Cauchy-Binet formula (12.1.71) 
in the last equality. 



Yl det(afc,iJ,<;,^^<^det(6fc,iJi 

jl>...>jm>l 

= det I ^aubki | /©. 

V i / l<k,i<m 

In particular for m = 1, 
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Hence 



det ( ^ akihi j 

V i ' l<k,e<m 



l<k/<m ' 

ending the proof of Lemma 12.61 ■ 

2.3 Probability on partitions expressed as Fredholm 
determinants 

Remember the definition fl2.0.3p of the Fredholm determinant of a kernel 
K{i,j), with i,j E Z and restricted to [n, oc). This statement has appeared 
in the literature in some form (see Case-Geronimo [28]) and a more analytic 
formulation has appeared in Basor-Widom [22]. The proof given here is an 
"integrable" one, due to Borodin-Okounkov |26j . 

Proposition 2.7 

P(A with Ai < n) = det ( / - K{k,- 



[n,n+l,...] 



where K{k, i) is a kernel 



K{k, 



2Tii ) J\^\=p<i J\z\=p-^>i z^^^w ^ z-w 

1 f iV IS dzdw z±V{z)-w£V{w) ^y^^^_y^^^ 



k — i \ 2-Ki J J J i»i=p<i z^^'^w'^ z — w 

z|=p-l>l 



for k ^ £, 

(2.3.1) 



with 



V{z) = -J2{t 



^ Z - SiZ }. 



The proof of Proposition 12.71 hinges on the following Lemma, due to 
Borodin-Okounkov ( [731 dl ES] ) : 
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Lemma 2.8 IfX = {xi, . . . , x^} C Z and S{\) := {Ai — 1, A2— 2, A3— 3, . 

P (A I S{X) D X) = i Yl ^A(t)sA(s) = det {K{x„ x,)), 



Xsuch that 
S(X)DX 



Proof of Proposition \2l[ Setting Ak := {X \ k e S{X)} with ^(A) := {Ai 



1, A2 — 2, A3 — 3, . . .}, one computes (the Xi below are integers) 

P(A with Ai < n) 
= P{X with all Ai < n) 
= P(A with5(A)n{n,n + l,n + 2, ...} = 0) 
= 1 — P{X with S'(A) contains some k for k > n) 

= 1- P{AnUAn+lU ...) 

n<i n<i<j n<i<j<k 

using Poincare's formula 

00 

= J2 ^(A with {a;i,...,x„}c 5(A)) 

m=0 n<xi<...<Xm 
00 

= ^(-l)'" J2 det{K{xi,x,))i<,,,<m 

m=0 n<xi<...<Xm 

= det (l~K{i,j] 



[n,n+l,...] 

from which Proposition 12.71 follows. 

Proof of Lemma \2. g[ " Remembering the probability measure introduced 
fl^XTD . we have that 

P(A|5(A)DX) = i ^A(t)3A(3). 

Asuch that 
S{A)DX 

Next, from the wedging-contracting operation 

^xC/a = /a, if a; G 

= 0, ifx^^A, 
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and using both relations fl2.2.7l) . one first computes: 

xex 

y xex A 

A \ x£X 

X 

S(X)DX 



^ fiCX 
S(X)DX 



= 1] SA(t)EsAv(s) (//^/a) = SA(t)SA(s). 

^ UCA such that 

S(X)DX S(X)DX 

Using this fact, one further computes (in the exponentials below the summa- 
tion E stands for Sf^) 



li that 
)DX 



Z 

X such that 
.S(X)DX 



z , 

using V'x,^^^. = -f'x^i^x, for i ^ j 

using (12:2:81) 

using Cauchy's identity f l2.2.4p 
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where 

Then, using Lemma [2.61 the expression above equals 

(^,„...^,,^:^...^:.^/0,/0> = det((^,,^:/0,/0))^^^^^^^ 

= det(ir(xfc,x^))i<^^^<„ 

where 
Using 

<^i^;/0,/0> = 1 if^=j<o 

= otherwise 

and setting 

V{z) = -J2iUz'' - Siz'), (2.3.2) 

i>l 

the generating function of the K{k, i) takes on the following simple form: 

</e-E*»A'e^^»^">(^)^*(u7)e-^^»^""eS*'^70,/0) 

. e- ^ *'^'e^ ^»^'>* («;)e- ^ ^«^""e^ *'^70, /0) 

eE^.^^e-^*>^"'e-^^''"'e^*>'"'' (V'(^)V'*(«^)/0, /0) using fl2X6D 
e^(^)-^('")(V^(^)^*(^)/0,/0) 
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,V(z)-V(w) 



w 



w 



- 1 + - + 



,V{z)-V{w) 



zJ 1-^ 

z 

,V(z)-V(w) ^ 



for \w\ < \z\ 



z — w 



and so 



K{k, 



1 

2Td 



\w\=p<l J\z\=p-'^>1 



1^1 z^~^^w~^ z — w 



ending the proof of Proposition 12.81 ■ 
Formula (12.0.51) in the remark is obtained by setting z ^ tz, w ^ tw in 

(TO, 



and taking 



d 



1 y 1 



- 



2m J t^- 



dz dw e^^*^)-^^*"') 



dt 



of both sides, using the t-independence of the left hand 



t=i 



side, yielding (I2.0.5P . 



2.4 Probability on partitions expressed as U(n) inte- 
grals 

This section deals with formula (I2.0.4p in Theorem 12. ![ stating that P(A [ Ai < 
n) can be expressed as a Unitary matrix integral. First we need a Lemma, 
whose proof can be found in |27j : 

Lemma 2.9 If f is a symmetric function of the eigenvalues Ui, . . . ,Un of 

the elements in U{n), then 

[ f = -J |A(^)|VK,...,n„)n-^. 
Ju{n) n\ J(5i)„ J-A 2muj 

Proof: Set G := U (n) and T := {diag(Mi, . . . , with Uk = e^^^ E S^}, and 
let t and q denote the Lie algebras corresponding to T and G. An element 
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in the quotient G/T can be identified with gT, because g' e gT imphes 
g~^g' e T, and thus there is a natural map 



Tx{G/T)^G:{t,gT)^gtg-\ 

Note that the Jacobian J of this map (with respect to the invariant measures 
on T, G/T and G) only depends out E i, because of invariance of the measure 
under conjugation. For a function / as above 

/ / = i / f{gtg-')J{t)dtd{gT) = X^M^ / f{t)J{t)dt. 
Jg n\ Jtx{g/t) n] Jt 

Denote the tangent space to G/T at its base points by t"*-, the orthogonal 
complement of t in = t © t"*". Consider infinitesimal changes t{l + e^) oit 
with ^ e t; also, an infinitesimal change of 1 e G, namely 1 i— > 1 + £77, with 
77 e t^. Then 



gtg-^-t = {l + er])t{l + e^){l-er] + 0{e''))-t 
= e{ti + r]t-tri) + 0{e^) 
=: £p + 0(£') 

and so 

Thus the Jacobian of this map is given by 

J(t) = det(yl(r^) -/), 

where A{t) denotes the adjoint action by means of the diagonal matrix 
denoted by t = diag(Mi, . . . , with \ui\ = 1. Let Ejk be the matrix with 1 
at the {j, k)-th entry and everywhere else. Then 

{A{t-') - I)Ejk = {ujW - l)Ejk 

and thus the matrix A{t^^) — / has n{n — 1) eigenvalues {uJ^Uk — 1), with 
i ^ j,k < n and j ^ k. Therefore 

det(^(t-^) - /) = nK"'«fc - 1) = n 1^^- - '^kW 

j^k j<k 

using the fact that uj^ — Uj. m 
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Proposition 2.10 Given the "probability measure" 

P{X) = Z-hx{t)sx{s), z = e^-i**'^% 
the following holds 

^ ' Juin) 

Proof: Using Lemma 12.91 and the following matrix identity, 

det (a,,<,(j) &j,aa))i<,j<„ = det (a»fc)i<i,fc<„ det (feife)i<,,fc<„ , 

one computes 



U{n) 

2iTizk J 

n / 



= I ydetfz^7\rrO T\{e-T.rit.^l^^^-C)A^ 

= nldetf / /-n.^-j:T{Uz^+s,z-^)Al\ =n!detm„(t,s) 
ending the proof of Proposition I2.10[ 



l<i,m<n 
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3 Examples 



3.1 Plancherel measure and Gessel's theorem 

The point of this section will be to restrict the probability 

P.,yW ■■= |s,(x)s,(|/) = |s,(t)s,(s) (3.1.1) 

considered in section 2 (formula fl2.0.1l) ) to the locus Ci, defined for real 
^ > (expressed in {x,y) and (s,t) coordinates), 

£i = {{x,y) such that ^xf' = = (5fciA/^} 

i>l i>l 

= {all Sk = = 0, except ti = Si = y^}. (3.1.2) 

The reader is referred back to subsection 11.5.11 for a number of basic formulae. 

Theorem 3.1 fj^ [13 W For the permutation group Sk, the generating 
function for distribution of the length of the longest increasing subsequence 

P\L{TXk) < n) 

is given by 



e-iy^lp\L{'Kk)<n) = e-«detf/ ^z^-^e^(^+^")^ 

= e-^det (^Jj_£(2v/^)' 



Ju{n) 



[n,n+l,...] 



where Ji{z) is the Bessel function and 

(4(2v/e)^m(2v/e) - J.+i(2v/e)^K2Ve)) , 

= J2Jk+n{2^/l)Je+n{2^/l). (3.1.3) 

n=l 
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Proof: For an arbitrary partition A G Y, and using formula fll.2.71) . the 
restriction of P^^y to the locus £i, as in (13.1.21) . reads as follows: 



P«(A) := P.,,(A) 



sxit)sx{s] 



^ |A|!^ |A|! 

-,^(/^ 
' |A|! |A|! 

= forn=|A|, 

where P(")(A) can be recognized as Plancherel measure on partitions in 
as defined in section 11.5.11 

P"(A) = A e Y„. 



n! 



It is clear that 



A G Y, 



|A|! V |A|! 

is a genuine probability (> 0), called Poissonized Plancherel measure. 
We now compute 



P«(A with Ai < 



n] 



AeY 

At <n 



^1^1 (/ 



A\2 



|A|! |A|! 



k\ ^ lAI! 

fe=0 |A|=fe I I 



fc=0 



A;! 
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and thus using Theorem I2.l|p| . 
P^(A with Ai < = Z-Met 



e"^ det 



dz 
51 2TTiz 

dz 



_^k-e^-j:^{tizi+s,z-i) 



l<k,e<n 



^k-l^-y^{z+z-^) 



51 2'Kiz 



l<kl<n 



U{n) 



e-«det (Jk-i{2^/^) 



l<k/<n 

where we used the fact that for a Toephtz matrix 



det(a^ ^Ck-e)i<k,e<n = det(cA:_^)i<fe,^<„, a 7^ 0. 



It also equals 



P«(A with Ai <n) = det ) 

V [n,n+l,...]J 

where K{i,j) is given by formula (12.0.51) . where V{z) = \^{z — z^^). Since 

^Tznz)-w£v{w) ^ f 1 _ ±y 

z — w \ wz J ' 



^The Bessel function J„(u) is defined by 

00 



and thus 



^(zz)"J„(2V^). 
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one checks 



{k-i)K{k,i) = (h (h dzdw' 



\z\=ci \ 
\w\=C2 



oo 

= (fc-£)^J,+42v^)J,+„(2v^). (3.1. 

n=l 

The last equahty follows from the recurrence relation of Bessel functions 

Jfe+i(2z) = -.h{2z)-.h-i{2z). 

z 

Indeed, subtracting the two expressions 

J,+i(2^) J,.(2^) = -Jfc(2^) J,(2^) - J,_i(2z) Jfc(2z) 
Jk+ii2z)Jei2z) = ^Jk{2z)M2z) - Jk-i{2z)J,{2z)., 

one finds 

(£ - k)Jk{2z)Je{2z) = z{Ju{2z)J,+i{2z) - Jk+i{2z)J,{2z)) 

-z{Jk-i{2z)Je{2z) - Jk{2z)J^^i{2z)) 

implying 

oo 

z{Jk{2z)J,+i{2z) - Jk+i{2z)J,{2z)) = {k-e)J2 Jk+n{2z)Ji+n{2z), 



n=l 



thus proving fl3.1.4p . 
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Remark: Incidentally, the fact that is a probability shows that Plancherel 
measure itself is a probability; indeed, for all ^, 

n=0 ■ AeY n=0 ' |A|=r!, ' n=0 ' AeY„ 

comparing the extremities leads to J^xeY^ P^'W — ^■ 
3.2 Probability on random words 

Here also, restrict the probability (I3.1.ip to the locus below, for ^ > and 

p G Z,>i, 



<^Xi = 6ki^, yi = ... = yp = p, all other yj = 



{tk = Ski^, ksk = pP''} 



Recall from section 1.5.4 the probability P'^'^ on partitions, induced from 
uniform probability P'^'^ on words. This was studied by Tracy- Widom [86] 
and Borodin-Okounkov |26j; see also Johansson |58j and Adler-van Moerbeke 

m 

Theorem 3.2 For the set of words S^, the generating function for the distri- 
bution of the length of the longest increasing subsequence, is given by (setting 
(3=1) 



-Ky^MLpfc,P(^(^^)<^) = e-f«detf / 

^' \JS' 2vriz / 

= e-P« I e^^'^ det{I + M)PdM 

Ju(n) 



det(l-K{j,k) 



(n,n + l,...) 
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with K{j, k) a Christoffel-Darboux kernel of Charlier polynomials: 
ij-k)Kij,k) 

dz dw \ p -j-j ^ 



kl=ci<j^ 
\w\=C2 

(1 - ^z)-Pe-^''' (1 - ^w)Pe^'"'' 



(P)j+ie^ 



(j + 1)! A;! 

{a)j := a{a + l)...{a + j - 1), (a)o = 1 
and iFi(a,c;x) is the confluent hypergeometric function; 



(3.2.1) 



2vr« Jj,|=,^<^^ m! 
— / {l-^wye^'"~'w'^-^dw = ^iFi(-g,m + l;e^). 

These functions are related to Charlier polynomial^ . 

Charlier polynomials P{k\a), with k G Z>o, are discrete orthonormal polynomials 
defined by the orthonormality condition 

oo ^ 

^P„(fc;a)Pm(fc;a)wa(fc) = for it;„(fc) = e"" — , 

with generating function 

oo 

a"/^—Pn{k; a)w" = £"""'(1 + w)^ 

n=0 
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Proof: This proof will be very sketchy, as more details will be given for 
percolation case in the next section. One now computes 



P«'P(A with Xi<n) = ^ e 



A 

Ai <n 



C2 



e-^«^5^s,(e,0,...)sA(p/3 



A 

Ai <ri 



E 

A-] <ri 



2 ' 3 

SA(l,0,...)S A(p,|,f,...) 

|A|! p 



|A| 



k\ ^ — ' p 

|A| = fc 



A;=0 



In applying Theorem I2.H one needs to compute 



/:2 



and 

Therefore 

P«'P(Ai < 



'C2 



£2 



Z-^ det 
e-«P^ det f 



dz 
51 27rzz 



l<kl<n 



C2 



\Jsi 2niz y i<ki<n 



-^'f^ det 



51 



dz 
2Tiiz 



l<k/<n 

the change of variable z 



U{n) 



e^^'^^ det{l + (3 MY dM. 
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Then, one computes 

z — w {1 — (3z){l — [3w) zw 

and 

leading after an appropriate rescaling, using f3 to formula (13.2. II) . From 
footnote [ini the confluent hypergeometric functions turn out to be Charlier 
polynomials. ■ 



3.3 Percolation 

Considering now the locus 

one is led to the probability appearing in the generalized permutations and 
percolations (sections 1.5.4 and 1.5.5), namely 

p{L{M)<i)= J2 (i-eni'isA(g,|,...)sA(p,|,...)- 

A 

\i < e 

We now state: (see [26l [56] ) 
Theorem 3.3 Assuming q > p, we have 
P{L{M) < e) 



Ju(e) 



det {I-K{i,3)) 
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where (C is a constant depending onp,q,C,) 

^ ^ Mpjp + ^)Mp.i{p + j) - Mp_i(p + z)Mp(p + 2) 

where 

Mp{k) := Mp{k;q-p+l,^) with k e Z>o. 
are Meixner polynomial¥^. 



Proof: Using the restriction to locus £3, one finds 
and 



= (1 - v^^)^(i - v^^-'r 



(1 - ^Tiz-y^i - 

V{z) = glog(l - ^-^) -plog(l - z). 



'^3 £3 

and thus 



Hence 



z-w ^i-^z){l-^w) zw{l-^z-^){l-Vlw-^) 

and so 



1 ^fn^)-<^H 



+ I. 2; — W 



-q 



^^Meixner polynomials are discrete orthogonal polynomials defined on Z>o by the or- 
thogonality relations 

and with generating function 

1 - 1) ^l - z)-^-^' = E /3, C). (3.3.1) 
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Expanding in Laurent series, one finds 
{l-r]z)-^{l-7]z-Y 



J2 + ^y^^ 



J=0 



^ + 1) ■ ■ ■ (/3 + 2 - 1) {-f3'){-[3' + 1) . . . {-(3' + J-1) 



i,j>0 



^ (?+m)! 7! ^ 

The coefficient of z"^ for m > reads 



2j+m 



00 



^ + 1) ■ ■ ■ (/3 + J + m - . . . i-P' + J - 1) 7]^^+"^ 

^ (j + m)! j! 

ml 

which is Gauss' Hypergeometric functior@ and thus 

One does a similar computation for the other piece in the kernel, where one 
computes the coefficient of z~'^. Using this fact, using a standard linear 
transformation formulaic for the hypergeometric function, and using a poly- 
nomial property for hypergeometric functionJ^ and the formula for < p < q 



1? 



r(l-g) ^ ' (q-p-iy. 



2Fi{a,b;c; z) :— ("-UW^ Notice that when a — — m < is an integer, then 

2^1(0, 5; c; z) is a polynomial of degree m. 

^•^In general one has 2^^1(0, c; z) = (1 — z)^'' 2Fi{b, c — a;c; jzij)- 

^"^When p is a positive integer (< q) and k an integer, then 2Fi{—p,—q;k + 

l;z) is a polynomial of degree p in z, which satisfies 2Fi{~p, --q; k + l]z) = 

r(fc+l)r(l-g+p) / F ( n k n-n -n 4- 1 • M 
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one finds Meixner polynomials (assume q > p > 0, x>0 and A; > 0) (see 
Nifikorov-Suslov-Uvarov [71J and Koekoek-Swartouw [65j): 

1 / 



■'\V^\<\z\=c^ 



(p+1) 



2Fi{p+l + k,-q;k + l;0 
e^'il-O' 2F,(^-p,-q;k + l;^^ 



{p + l)kT{l+p-q) e^^l^ , 

(fc + p)!r(i-g) (1-0^'^ - ^ - ^ + 1' 

r(l+p-g) ^p+'^z^ 



p\T(l-q) {l-i)P-'' 



—Mp{p + k;q-p+l;0 



P J (1-0 



P-9 



Mp(p + /c;g-p+l;0, 



where 



Mp{x;P,0= 2F,i-p,-x;P;^) = (^1^^ J^x^ + . . . , x G Z>o 

=: QpX^ + . . . , with P = q — p + 1 

are Meixner polynomials in x, satisfying the following orthogonality proper- 
ties: 

OD , 

w{x)Mp{x](3,C)Mm{x](3,C) = ^p^-^_^^^^^^ Spm =■■ hp6pm, 
with weight 

wix) := %^ = il±lf^e^ 
a;! 
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Using these facts, one computes 



p 



(27ri)2 J J ivei<Ni=c 



dz dw 



w\=C2<\ 



(1 - z)-^-\i - z-y (1 - ^r'(i - ^-') 



n-5 



-9 



(1 - zr^ii - VI z-'y-' (1 - wy{i - VI w-') 



l\-q-l 



yk+2 



w 



-e+1 



= -^-^VHp + k)w{p + e) {Mp{p + k)Mp_r{p + £)- Mp_,{p + k)Mp{p + £)) , 
where 

Mp{p + k) := Mp{p + k-q-p + l,0- 
We also have, using Propositions 2.2 and 2.10 restricted to the locus £3, 

P(L(M) < t) 



- E (i-eni^'sA(?,|,---)sA(p,|,...) 



X 



kh = qe'^ 
ksk = p^l"^ 



dz \ 



= (l-erdet(^^ z'-'-'e-^Titjz^+sjz 

= (1 - ^y^det(^^z--'^'ii-V^zy{i-V^z-'y^^ 
= (izi)!! f \Ae{z)\'fl{i+V~^z,y{i+V~rz,r 

= / det(l + ^/^My det{l + y/^MydM, 

Ju(e) 



ktk = qe^^ 
ksk^pi^'"^ 



l<a,a'<€ 



dZj 



27rizj 



'u{(.) 

thus ending the proof of Theorem 3.3. 
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4 Limit theorems 



4.1 Limit for Plancherel measure 

Stanislaw Ulam [88] raised in 1961 the question: how do you compute the 
probabihty P*^"^(L(7r„) < k) that the length L := L^"'^ of the longest increas- 
ing sequence in a random permutation is smaller than k. What happens for 
very large permutations, i.e., when n — > cxo? By Monte-Carlo simulations, 
Ulam conjectured that 

c := lim ^ (4.1.1) 



exists (i?*^"^ denote the expectation with respect to P^^^). A much older 
argument of Erdos and Szekeres [36] implied that E^'^\L) > \\/n — 1 and 
so c > \. Numerical computation by Baer and Brock [17] suggested c = 2. 
Hammersley [50j showed the existence of the limit (14.1.11) : in 1977, Logan and 
Shepp [nZ] proved c > 2 and, at the same time, Vershik and Kerov [22] showed 
c = 2. More recently other proofs have appeared by Aldous and Diaconis 
[15], Seppalainen [78] and Johansson [57]. Meanwhile, Gessel [15] found a 
generating function for the probability (with respect to n) and connected this 
problem with Toeplitz determinants. Monte-Carlo simulations by Odlyzko 
and Rains [72] suggested 

VarL(") , E(")(L(")) -2Vn 
Co ~ lim — — ci = hm ^ , 

with Co ~ 0.819 and Ci 1.758. 

In this section, we explain some of the ideas underlying this problem. It 
is convenient to write a partition A G Y in (m, t>)-coordinates, as shown in 
the picture below. 
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Figure 4.1 

Remember from (11 .2. 51) the hook length h^, defined as 

/i^ = JJ h^j, with h^j = hook length = Aj + Aj - z - j - 1, 

where Aj is the length of the j^^ column. Also remember from ( 11.2.60 the 
formula expressed in terms of h'^. Then one has the following Theorem, 
due to Vershik and Kerov [9lj; see also Logan and Shepp [67j. A sketchy 
outline of the proof will be given here. 

Consider Plancherel measure (see section [3?T!) . which using formula (ll.2.6p . 
can be written in terms of the hook length, 

and define the function 

{|(m arcsin| + V4 — m^), for \u\ < 2 
(4.1.2) 
for \u\ > 2. 
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Theorem 4.1 ( Vershik-Kerov-Logan-Shepp) l91\ Upon expressing A in 
{u,v)- coordinates, define the subset of partitions, for any e > 0, 



YJe) :-- 



A G Y I sup 



n 



X{uy/n) — VL{u) 



< e 



Then, for large n, Plancherel measure concentrates on a Young diagram 
whose boundary has shape Q{u); i.e., 

lim P(")(Y„(£)) = 1. 

n— >oo 

Moreover, for the length Ai of the first row, one has 

Ai 



lim P(") 



2v^ 



1 



< £ 



and thus for uniform measure on permutations (remembering L^iTn) = the 
length of the longest increasing sequence in iTn), one has 



lim P(") 



L(tt„ 



2v^ 



< e 



1. 



Brief outline of the Proof: In a first step, the following expression will be 
estimated, using Stirling's formula^ 



— logP„(A) 

n 

1 , n\ 

= log ■ 



-log n 

n J. J. 

(i,i)eA 

n ^i- 



— log n! 

n 



(i,i)GA 



logn log v27m + 1 + 0(1), by Stirling 

n 

l + -(log TT hf--nlogn}'^^) log V27rn + 0(1), using lAI 

(ij)6A 

1 + 2 Y: -loF^^' — 



n, 



+0 



(Riemann sum) 



15 



log nl — n log n + log \/2TTn — n 



for n y oo. 
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-^1 + 2 dx dy\og{F-^{x) -y + F^\y) - x), 

J {(x,y),x,y>0,y<F^(x)} 

assuming the partition X = ^ tends to a continuous curve y = Fj^{x) in 
(x, y) coordinates. Then the Riemann sum above tends to the expression 
above for n — > oo; note 1/n is the area of a square in the Young diagram, 
after rescahng by and thus turns into dxdy in the hmit. Clearly, the 

expression Fj^{x) — y + F^^ (y) — a; is the hook length of the continuous curve, 
with respect to the point {x,y). In the (w, vj-coordinates, this hook length 
is particularly simple: 

F^{x) - y + F^^{y) - x = \Pl{u2 - with u\ < U2, 

where (i^i, = (^)-^a(^)) {u2,'ijj{u2)) — {F^^{y),y) and where v — 

ip{u) denotes the curve y = Fj^{x) in (m, t;)-coordinates. 

Keeping in mind Figure 4.1, consider two points {ui, iIj{ui)) and {u2, 'ip{u2)) 
on the curve, with ui < U2. The point {x,y) such that the hook (emanating 
from (x, y) parallel to the (x, |/)-coordinates) intersects the curve at the two 
points (rti, -0(1*1)) and (1*2, ^'('^2)), is given by a point on the line emanating 
from (wi,'0(iii)) in the (1, — l)-direction; to be precise, the point 

1 

{x, y) = {ui, ip{ui)) + -{u2- ^(^2) - Ml + i'{ui)){l, -1). 

So the surface element dxdy is transformed into the surface clement 
duidu2 by means of the Jacobian, dxdy — 1(1 +'0'(iii))(l — ip'i.U2)duidu2, 
and thus, further replacing ui 1— > V^Ui and U2 ^ V2u2, one is led to 

-ilogP„(A) 
n 

~ 1 + 2// dx dy\og{F{x) + F+{y)-x-y) 

J J {{x,y),x,y>0,y<F{x)} 



= 1 + 



/ / duidu2{l + '0'K))(1 - ^''(^^2)) log{V2(u2 - Ui)) 

J Jui<U2 

— - 11 log{V2\u2 — ui\)f{ui)f'{u2)duidu2 + 2 f f{u) a,TCCOsh.\^\du 
2JJm2 J\u\>2 2 

1 rr r fM^fM)\u,du2 + 2 / /(«) arccoshl^^i. 

2 77r2 V U1-U2 ) J\u\>2 2 
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upon setting ^ip^u) = Q{u) + f{u), where the function Q{u) is defined in 
fl4.1.2l) . The last identity is obtained by using Plancherel's formula of Fourier 
analysis, 

9i{v)g2{v)dv = — / gi{v)g2{v)dv 



2vr 

applied to the two functions 

r p -F f \ 

gi{v) := / dui\og{V2\v - ui\)f'{ui) = / — — -dui 

Jr jRy/2\v-Ui\ 

and g2iv) := f{y). 

This shows the expression —y^ logP„(A) above is minimal (and = 0), when 
f{u) = 0; i.e. when the curve ipi"^) = ^{u) and otherwise the integrals above 
> 0. So, when the integral equals e > 0, the expression — ^ logP„(A) ^ e > 
and thus P„(A) ~ e"*^", which tends to for n y oo. Only, when £ = 0, is 
there a chance that P„(A) = 1; this happens only when ip = Q. ■ 



4.2 Limit theorem for longest increasing sequences 



In section 13.11 it was shown that a generating function for the probability of 
the length L{'Kk) of the longest increasing sequence in a random permutation 
is given in terms of a Bessel kernel: 



e-^Y.-^P\L{Tiu)<n) = det{l-K{j,i)) 

k=0 



[n,n+l,...] 



where 

oo 

K{k,e) = Jk+U^V^)Jumi2^). (4.2.1) 

m=l 

In the statement below A{x) is the classical Airy function 

A[x) = — cos \- xu ] du; 

Jo V 3 / 

this function is well known to satisfy the ordinary differential equation A"{x) = 
xA{x) and to behave asymptotically as 

_ 2 3,3/2 

^(^) ■■= 1/4 + as X ^ oo. 
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Theorem 4.2 (J. Baik, P. Deift and K. Johansson /7P|/) The distribution of 
the length L{TTn) of the longest increasing sequence in a random permutation 
behaves as 

lim P(") (L(7r„) < 2n'/^ + xn^/^) = det (J - Ax[.,oo)) , 



where 

A{x,y) :-- 



, ^ X A(x)A'(y) - A'(x)A(y) 

duA{x + u)A{y + u) = ^ ' K_±^_ ^4 2. 2) 

x-y 



The proof of this Theorem, presented here, is due to A. Borodin, A. Okounkov 
and G. Olshanski |25]; see also [20]. Before giving this proof, the following 
estimates on Bessel functions are needed: 

Lemma 4.3 (i) The following holds for r ^ 00, 

k'/'^2.w/3(2r)-A(a;)| =0(r-i/3), 

uniformly in x, when x G compact K cM.. 

(ii) For any 6 > 0, there exists M > such that for x,y > M and large 
enough r, 



s=l 



< 



Lemma 4.4 (de-Poissonization lemma) (Johansson 151 j) Given 1 > Fq > 
Fi>F2> ... > 0, and 

00 ^i. 



there exists C > and k^, such that 



F{k + A^k\ogk) - ^<Fk<F{k- A^k\ogk) + ^, for all k > ko. 
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Sketch of proof of Theorem 4-2: Putting the following scaling in the Bessel 
kernel K{k,£), as in fl4.2.ip . one obtains, setting r := 

N oo 
k=l k=N+l 



N oo 

(E+ E) 

h=\ k=N+l 

^-1/3 [ri/V2,+(,+fc,-i/3),i/3(2r)] [r^/V2,+(j,+fc,-i/3),i/3(2r)] (4.2.3) 



Fix S > and pick M as in Lemmall^(ii). Define N := [(M - m + l)r^/^] = 

O(r^), where m is picked such that x,y >m (which is possible, since x and 
y belong to a compact set). Then 



k=N+l 



J2r+ri/3(a.+fcr-l/3)(2r) J2r+rV3(y+fcr-l/3)(2r 
oo 

J2r+rl/3(a;+Af_m+l)+s(2'")'^2r+rl/3(y+Af_m+l)+s(2'") 



(4.2.4) 

the latter inequality holds, since 

X + {M - m + 1) = M + {x - m) + 1 > M. 
On the other hand, 

|r^/V2,+,i/3(,+fc,-i/3)(2r) - A{x + kr-'/^)\ = 0{r-'/^) 

uniformly for x G compact K G and all k such that 1 < k < N = 
[(M — m + l)r^/^] . Indeed, for such k, 



m 



<x<x + kr < X + M - m + 1 = M + (x - m) + 1 



and thus, for such k, the x + kr ^/^'s belong to a compact set as well. Since 
the number of terms in the sum below is = [(M — m + l)r^/^]. Lemma H3] 
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(i) implies 



N 

|Xl(^^'^2r+xri/3+fc(2r)r3j2^^^^i/3+fe(2r) 

k=l 



-A{x + kr-^/^)A{y + kr^/^)) =0(1). (4.2.5) 



But the Riemann sum tends to an integral 

[(A/-m+l)ri/3] 

^-1/3 



V A{x + kr-^/^)A{y + kr~^/'^) ^ / A{x + t)A{y + t)dt. 

k=i Jo 

(4.2.6) 

Hence, combining estimates (14.2.41 14.2.51 I4.2.6P and multiplying with r^/^ 
leads to 



(2r)J2r+„ri/3+fc(2r) - 



M-m+1 



A{x+t)A{y+t)dt 



< 6+o{l). 



k=l 

for r oo. Finally, letting 5^0 and M — > cxd leads to the result. Thus the 
expression (14.2.31) tends to the Airy kernel 

POO 

A{x,y) := / duA{x + u)A{y + u). 
Jo 

for r = — > CX3. Hence 

lim e-« V (L(7r„) < 2^^/^ + x^^'') 
f^oo ■^^ — ' n! 

n=0 



lim det / - K{i, t) 

^— >oo 



[fc,fc+l, 



A; = 2^1/2 + x^i/e 



oo „ k 

= 1 + V(-1)M det(A(z„z,)) TTdz, 

= det (/ - Ax[x,oo)) • 

Finally, one uses Johansson's de-Poissonization Lemma 14.41 From Corol- 
lary [L5l Plancherel measure Pn(Ai < xi,...,Afc < Xk) decreases, when n 
increases, which is required by Lemma 14. 4[ It thus follows that 

lim P (L(7r„) < 2n^/^ + xt/"") = det (/ - Ax[.,oo)) , 

n— ►cxD 

ending the proof of Theorem 14. 2[ ■ 
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4.3 Limit theorem for the geometrically distributed 
percolation model, when one side of the matrix 
tends to oo 

Consider Johansson's percolation model in section ll.7.2[ but with p and q 
interchanged. Consider the ensemble 

Mat^^'P) = {qxp matrices M with entries = 0, 1, 2, . . .} 

with independent and geometrically distributed entries, for fixed < < 1, 

P{M,, = k) = {l-i)i\ A; = 0,l,2,.... 



{over right /down 
y Mij, paths starting from 
' ^ / \ / \ 

entry (1, 1) to 

has the following distribution, assuming q < p, 

p{L{M)<i) = z-i \ihu...,Krf[^J^^±^p)^e^ 

hem i=i *■ 

max(/i^ ) <-^+t? — 1 

where 

9-1 

Zp,, = ^"^(1 - O-'W U^Kp - q + j)! (4.3.1) 

i=o 

Assuming that the number of columns p of the g x p random M matrix 
above gets very large, as above, the maximal right/lower path starting from 
(1,1) to (g,p) consists, roughly speaking, of many horizontal stretches and q 
small downward jumps. The Mij have the geometric distribution, with mean 
and standard deviation 



k=0 1 ^ 

a^(M,,) = E(M^.)-(^(M.,))^ = 5^PP(M,,=fc)-(^y 



e(e + i) f e 



1-0' vi-e; (1-0 



2 
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So, in the average, 

L{M) ~ pE{Mij) = for p ^ oo 

with 



' (1-0'' 
Therefore, it seems natural to consider the variables 

_ L(M) - _ Ai - 

and, since A = (Ai, . . . , A,), with q finite, all Aj should be on the same footing. 
So, remembering from the proof of Theorem 11.91 that hi in (14.3. ip and Aj are 
related by 

hi = q + Xi- 



we set 



_ - Yz^p 



Theorem 4.5 (Johansson l5E]) The following limit holds: 



1 



lim P — — ^ < y 



Jm 1 



which coincides with the probability that a q x q matrix from the Gaussian 
Hermitian ensemble (GUE) has its spectrum less than y; see section V8.2.1\ 
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Proof: The main tool here is Stirhng's formula^. Taking into account hi = 
q + \i — i, we substitute 

Xi = Xi 

in (14.3.11) . The different pieces will now be computed for large p and fixed q, 

Q 

Y[{hi+p-q)\ 

1=1 

g 

1=1 



1=1 



2 log - 1 + log 1 + a;i ' 



^ f P , A 



1 — ^ V y p p 



C 2 IP I o M 

^27rp^^ "^/^ ^2(1 - V l-r KVp) 

QP , + 1) ^ 1 P 

e V ^ ^1 / M 1 + 



Stirling's formula: 

n! = %/2^ e"('°s"~i) fi + O fi^ ^ , for 
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upon expanding the logarithm in powers of 1/ ^/p. Similarly 

i=l 

i=l 

■ mm 



Xj ^ (g- 0(1-0 



i + o{- 

p 



.e 



.e 



.e 



, Vp^ sr^ . 
1 / Xi + 



log 



1-e 



Also 



logC 
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Besides 

9-1 9-1 / _ A / /I \ \ 



j=0 i=0 

1-1 



(2.p)'/^e-" + — f 1 + O ^ ' 



and 



Ag{hi,...,hq) = Yi i^i-hj) 

l<i<j<q 



l<i<i<g 



1-e 



Remembering the relation hi = q + Xi — i, and using the estimates before in 
equahty = below, one computes 



hem 1 

max hi<.£.-\-q—l 



max 1 
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[ Aixyf[e-^"/^dx, 

J Xi<y ^ 

q ■ 



using Selberg's formula 



and noticing that in 

hi = q + Xi-z = YZ7^ + r^^^^ + q- i 
an increment of one unit in Aj implies an increment of Xi by 

Therefore, one has 
Finally, for p large, one has 



1 



hi<i + q-i = + + 9-1 < y- 

The connection with the spectrum of Gaussian Hermitian matrices will be 
discussed in section 18.2. 1[ This ends the proof of Theorem 14. 5[ ■ 



4.4 Limit theorem for the geometrically distributed 
percolation model, when both sides of the matrix 
tend to oo 

The model considered here is the percolation model as in section 14. 3[ whose 
probability can, by section 13.31 be written as a Fredholm determinant of a 
Meixner kernel; see Theorem 13.31 
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Theorem 4.6 (Johansson l5E]) Given 



i/ie following holds: 



lim 

p,g -> oo 
g/p = 7 > 1 fixed 



PP 



1/3 



which is the Tracy- Widom distribution (see section W^) . In other terms, the 
random variable L{Mg^p) behaves in distribution, like 



i 

oo 



oo 
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Sketch of Proof : Remember from Theorem I3.3[ the probabihty P{L{Mq p) < 
z) is given by a Fredholm determinant of a Christoffel-Darboux kernel com- 
posed of Meixner polynomials. Proving the statement of Theorem 14.61 amount s 
to proving the limit of this Meixner kernel with the scaling mentioned in the 
Theorem tends to the Airy kernel, i.e., 



lim pp^^^Kp{ap + 
wherj^ 

Kp{x,y) = - 



pp ' ^r],ap + pp '^TTj ) = Mrj.rj ) := ; 

7] — 1]' 



x + P'-l\ fy + p'-l 



X 



mpix, P', i)mp_i{y, p', ^ - mp{y, p', ^)mp.i{x, /?', 



x-y 



for the Meixner polynomials, which have the following integral representa- 
tion, a consequence of the the generating function ( 13.3.11) . (slightly rescaled) 



mp{x,P',0 =p\C 



— dz 



zP{l - zY+P' 2mz 
Thus, it will suffice to prove the following limit 



lim V ' 



{y - 0" 



dz 



where 



\z\=r<i zP{z - l)^+[7p]-p+i 2mz 



CA{r]), (4.4.1) 



x=f3p+pp^^'^ri 



1-e ' /? = «-7 + l= 



p \ I It 



i^ith dp ^ (i4)'tg-^l:i), and P' = q-p+l^p^-p+l. 



In view of the saddle point method, define the function Fp{z) such that 



2:P(^Z — l)^+[7P]-P+l 

Then one easily sees that 



F,{z) = F{z) + pp-^'S log ^ + (7 - ^) log(^ - 1), (4.4.2) 

z — I p 

where the p-independent function F{z) equals, 

F{z) = (3 log{z - - « log(^ - 1) - log (4-4.3) 

For the specific values above of a and (3, one checks the function F{z) has a 
critical point at := — a/C/T) i-e., 

F' (z,) = F" (z,) = 0, F'" (z,) 



e(v^+v^)(l + v^) 
and thus 

F{z) - F (z,) = kz- z^fF'" {z,) + 0{{z - z^Y). (4.4.4) 



Setting 

z = z,{l - p-^/lsC)e''^P~'^' (4.4.5) 

in Fp{z) as in fl4.4.2p . one first computes this substitution in F{z), taking 
into account (14.4.41) . 

F{z) = F(^J + !iz|£^(t + ,,)3ir'"(^j + 0(p-4/3) 

bp 

= F{z,) + ^{t + zsf + 0{p-'/'), (4.4.6) 
Sp 

upon picking -in the last equality- the constant C such that 

{-z,C)'F"' (zj = 2. 

Also, substituting the z of fl4.4.5p in the part of F{z) (see fl4.4.3p ) containing 

V, 
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9/-? , z — £ or) ^ Zc — £ pCzc{£ — I) iriit + is) ^, .in. 

Z - 1 p2/3 _ 1 _ ^) p 

(4.4.7) 

Thus, adding the two contributions (14.4.61) and (I4.4.7p . one finds 

pFp (z,(l - p-i/3^C)e**^^'"'') = pF,{z,)+'-(^it+tsf+3vit+ts)'^ +0 {p-'/'') 

(4.4.8) 

One then considers two contributions of the contour integral about the circle 
\z\ = r appearing in (14.4. ip . a first one along the arc (vr — 5p, vr + 6p), for 
6p tending to with p oo and a second one about the complement of 
(tt — 6p,7i + 6p). The latter tends to 0, whereas the former tends to the Airy 
function (keeping s fixed, and in particular = 1) 



/'OO 3 

2vr Joo 



upon noticing that dz/z = iCp^^^^dt under the change of variable z t, 
given in (14.4.51) . establishing limit (14.4. ip and finally the limit of the Meixner 
kernel and its Fredholm determinant. Further details of this proof can be 
found in Johansson [56]. The Fredholm determinant of the Airy kernel is 
precisely the Tracy- Widom distribution, as will be shown in section 19.21 ■ 



4.5 Limit theorem for the exponentially distributed 
percolation model, when both sides of the matrix 
tend to cxD 

Referring to the exponentially distributed percolation model, discussed in 
section [T. 101 we now state 

Theorem 4.7 (Johansson 156 j ) Given 



(1 + vir^' 
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the following limit holds: 



lim PI + - 

I (pl/2 + gl/2)(p-l/2 + ^-1/2)1/3 - 



g/p = 7 > 1 /ixed 

1™ — ^^^v^ — <?/ =^(w), 

p,q ^ oo \ HF / 

q/p = 7 > 1 /ixed 

which is again the Tracy-Widom distribution. Here again L[Mq^p) behaves, 
after some rescaling and in distribution, like the Tracy- Widom distribution, 
for large p and q such that q/p = J > I.' 

L{Mg^p) r^ap + pp^l^r 

Proof: In Theorem 11.101 one has shown that P{L{M) < t) equals the ratio 
of two integrals; this ratio will be shown in section [7| on random matrices 
(see Propositions 17.91 and 17.41) to equal a Fredholm determinant of a kernel 
corresponding to Laguerre polynomials: 

J[Q,oo]i py^) iii=i^j ^ 

where 



hp-1 X - y 

1/2 



in the formula above, the = "^^3;" + . . . = (—1)^ ((^^Lyi) -^p(^) 

are the normalized Laguerre polynomial^: 







18^1") (y) = E:=o(-1)"(:-:) S = ^/cS^'^^. where C is a circle about 
z = \. 
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Therefore, using the precise values of a and p above, 

P{^^^^^^jf^^y) = det(/-/C(e,r,)X[.oo)) 

with 

/C(e, T]) = bp^/^Kl^-^^P {ap + pp^'^i, ap + pp^'S) 
The result follows from an asymptotic formula for Laguerre polynomials 

lim K^^-^'^'P [ap + pp^'^i, ap + pp^'^rj) = A(^, r/). 



with A{x,y) the Airy kernel, as in fl4.2.2l) . namely 

Aix)A'iy) - A'ix)Aiy) 



A{x,y) :-- 



x-y 



The Fredholm determinant of the Airy kernel is the Tracy- Widom distribu- 
tion. This ends the proof of Theorem 14. 7[ ■ 



5 Orthogonal polynomials for a time depen- 
dent weight and the KP equation 

5.1 Orthogonal polynomials 

The inner-product with regard to the weight p{z) over M, assuming p{z) 
decays fast enough at the boundary of its supporO 

{f,g)= [ f{z)g{z)piz)dz, (5.1.1) 
Jr 

leads to a moment matrix 

= {^'ij)o<i,j<n = ( (^*) 2:"' ) )o< j j<n-l • (5.1.2) 

Since the pij depends on i + j only, this is a Hdnkel matrix, and thus sym- 
metric. This is tantamount to the relation 

Amoo = rriooA^ , 

^^In this section, the support of the weight p(z) can be the whole of R or any other 
interval. 
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where A denotes the semi-infinite shift matrix 



/ 1 ... \ 
1 ... 
1 ... 



V 



Define 

Tn := detm„. 

Consider the factorization of rrioo into a lower- times an upper-triangular 
matrix0 

moo = S-^S^-\ (5.1.3) 

with 

S = lower triangular with non-zero diagonal elements. 
For any z G C, define the semi-infinite column 



(5.1.4) 



and functions Pn{z) and qn{z) 



Pn{z) := {Sx{z)),, 

Quiz) := {S^-\{z-'))^. 



(5.1.5) 



(1) The Pn{z) are polynomials of degree n, orthonormal with regard to p{z), 
and qn{z) is the Stieltjes transform of pn{z), 

Jr z-u 

Indeed, 

{{Pk^Pi)) = / Sx{z){Sx{z)yp{z)dz = SmooS'^ = I. 

\ / o<k,e<oD Jjj 

^°This factorization is possible as long as t„ := detm„ ^ for all n>l. 
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Note that Sx{z){Sx{z))'^ is a semi-infinite matrix obtained by multiplying 
the semi-infinite column Sx{z) and row {Sx{z))~^ ■ The definition of 
together with the decomposition 5''''"^ = Srrioo, leads to 

Quiz) = {S^-\{z-')^ 

j>0 

j>0 o<e<n 



o<e<n j>0 
z — u 



(2) The orthonormal polynomials Pn have the following representation 



Pn{z) 



1 



■det 









rrin 




Z 


\ l^nfl ■ ■ ■ 







(5.1.6) 



As a consequence, the monic orthogonal polynomials Pn{z) are related topn{z) 
as follows: 



Pn{z) 



Tn+l 



-Pn{z). 



(5.1.7) 



Defining p'n{z) to be the polynomial on the right hand side of (15.1.61) . it 
suffices to show that for k < n, 



(p'n^p'k) = and {p^,pj = 1, 
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thus leading to p„ = 


Indeed, 






/ 




(K, z^) = ^ det 








\ /^n,0 • • • fJ'n,n-l 


{z'',z'') J 



0, for k < n, 



for k = n, 



and thus for k = n. 



from which (15.1.61) follows. Formula fl5.1.7p is a straightforward consequence. 

(3) The monic orthogonal polynomials pn and their Stieltjes transform have 
the following representation 

det {fiij-lfii,,+i)^^. .^^_^ 



Pn{z) = Z""- 
Pn{u)p{u) 



-du 



det (A''ji)o<i,j<n~i 

^det (/ii, + + jl^i,j+2 + . . •)o<i,i<n 



z ~ u 



det (/i. 



»J/0<i,j<n-l 



(5.1.8) 



Proof: Setting 



one computes 

det (/iy - 2;" Vij+i) 



0<i,i<n-l 

det (2;/iij - /^ij+i)o<ij<„_i 

det(2;/I([ - - /^L ^/^I-i - P'l) 



det I E E ^^^^ii:^, znU - III 



Z3 



by column operations 



det ( zfx^ - - zfil_^ - /I^ 
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— det 



1 , 

— det 



( zflo 



\ 







enlarging the matrix by one row and column 

z \ 



= TnPniz). 

Setting this time 



by adding a multiple of the last row to rows 1 to n 



one computes 



pn+l 



det(/i,, + ^ + ^ + ... 

(OO -Tf OO -Tf- OO 
OO -.T 



0<i,j<n 



det jl'^ , jlj , . . . , ill_i,Y^ 



Z3 



OO 



z^ det I • • • , XI ^ 
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det 



oo ^ 

J2 (^) Piu)du 

j=0 

oo J 

^— j up{u)du 

3=0 ^ 



oo 



oo ^ 
j=0 ^ 



5 A'-n-l 



1 \ 

U 



u 



» oo 



z — u 



du 



Remark: Representation ( IS.l.Sp for orthogonal polynomials p„ can also be de- 
duced from Heine's representation. However, representation (15.1.81) is much 
simpler. 

(4) The vectors p and q are eigenvectors of the tridiagonal symmetric matrix 

L := SAS-\ (5.1.9) 

Conjugating the shift matrix A by 5 yields a matrix 

L = SAS-^ 

= SAS-^S^-^S^ 

= SAniaoS^ , using (I5.1.3p . 

= S'mooA^S'^, using Amoo = mooA^, 

= S{S-^S^-'^)A'^S'^, using again fl5X3|) . 

= {SAS-^y = L^, 

which is symmetric and thus tridiagonal. Remembering x{^) = (I5 ^) ^"^^ •••)^) 
and the shift (At>)„ = Vn+i, we have 

^Xiz) = zxiz) and A'^ xi^"^) = ^xiz'^) - zci, with ci = (1, 0, 0, ...)'^. 
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Therefore, p{z) = Sx{z) and q{z) = S'^ ^xi^ ^) are eigenvectors, in the 
sense 

Lp = SAS^^Sx{z) = zSx{z) = zp 
L^q = S^-^A^S^S^-\{z-^) = zS^-\{z-^)-zS^-^ei = zq-zS^-'ei. 

Then, using L = , one is led to 

((L — zl)p)n = 0, for n > and ((L — zl)q)n = 0, for n > 1. 



(5) The off- diagonal elements of the symmetric tridiagonal matrix L are 
given by 

Ln-l,n=\h^. (5.1.10) 

Since {pn,Pn) = hn, one has Pniy) = -k=Pniy)- From the three step relation 
Lp{y) = yp{y), it follows that 

?/" + ...) = ?/p„_i(?/) = L„_i,„p„(?/) + (terms of degree < n - 1) 



Ln-l,n ( ^?/" + • • • 



leading to statement (IS.l.lOp . 

5.2 Time dependent Orthogonal polynomials and the 
KP equation 

Introduce now into the weight p{z) a dependence on parameters t = (ti, ^2, • • •): 
as follows 

ptiz):=piz)e^^'^^\ (5.2.1) 

Consider the moment matrix mn{t), as in (15.1.21) . but now dependent on t, 
and the factorization of moo into lower- times upper-triangular t-dependent 
matrices, as in (15.1.31) 

moo(t) = (/iy(t))o<*j<oo = S-\t)S^^\t). (5.2.2) 
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The Toda lattice mentioned in the Theorem below will require the following 
Lie algebra splitting 

gl{n)=s®b, (5.2.3) 

into skew- symmetric matrices and (lower) Borel matrices. 

Also, one needs in this section the Hirota symbol: given a polynomial 
p{ti, t2, ...) oi a. finite or infinite number of variables and functions /(ti, ^2, • • •) 
and g{ti,t2, . . .), also depending on a finite or infinite number of variables ti, 
define the symbol 

The reader is reminded of the elementary Schur polynomials e^i°*'^' : = 
^^>g Sj(t)z* and for later use, set for £ = 0, 1, 2, . . ., 

One also needs Taylor's formula for a C°°-function /: 

f(z + y)=ey^^fiz), (5.2.6) 

which is seen by expanding the exponential. The following Lemma will also 
be used later in the proof of the bilinear relations: 

Lemma 5.1 1^13^ If denotes the integral along a small circle about oc, 
the following identity holds (formal identity in terms of power series): 

f{u)g{u)du =^<f dzf{z) [ ^du, (5.2.7) 

-^^^ Joe Jr ^ ~ U 

for holomorphic f{z) = Yli>o ^i^^ '^^'^ 9i^)' latter assumed to have all its 
moments. 

Proof: For holomorphic functions / in C, 

/ dzfiz) [ ^du = Res,=oo(Va.^^)(- V^-^' / g{u)u^du) 

= / g{u)u^du 

g{u) ttiU^du 

■ j>0 

{f,9), (5.2.8) 
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ending the proof of Lemma 15.11 



The next Theorem shows that the determinant of the tzme-dependent 
moment matrices, satisfies the KP hierarchy, a non-hnear hierarchy, whereas 
in the next section, it will be shown that these same determinants satisfy 
Virasoro equations. These two features will play an important role in random 
matrix theory. Notice that this result is very robust: it can be generalized 
from orthogonal polynomials to multiple orthogonal polynomials, from the 
KP hierarchy to multi-component KP hierarchies; see [12j . 

Theorem 5.2 0, The determinants of the moment matrices, also repre- 
sentable as a multiple integral, 



Tn{t) 



I r 

:= detmnit) = - / Al{z)Y\p{zk)e^^-^'^''>'dzk (5.2.9) 



satisfy 

(i) Eigenvectors of L: The tridiagonal matrix L{t) admits two independent 
eigenvectors: 

{L(t)p(t; z))n = zpn{t;z),n>0 
{L{t)q{t]z))n = zqn{t]z), n>l. 

• Pn{t', z) are nth degree polynomials in z, depending ont E C°°, or- 
thonormal with respect to pt{z) (defined in Ii5.2.1\) ). and enjoying 
the following representations: (define x{z) := {1, z, z^ , ...)~^ ) 

p^{t; z) := {S{t)x{z))n = ^X-^/^ ^"^^~lf'^^ , ■■= 

Tn[t) Tn[t) 

(5.2.10) 

• qn{t, z), n > 0, are Stieltjes transforms of the polynomials Pn{t', z) 
and have the following r-function representations : 

qn{t-z):=z! ^^p,{u)du = {S^-\t)x{z-')) ^ 



u 

^~n.~l/2 ^n+l{t+[z-^]) 
Tnit) ■ 

(5.2.11) 

^A„(z) = ^l<^<J<„(^^-^i) 
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(ii) The standard Toda lattice , i.e., the symmetric tridiagonal matrix 



L{t) := S{t)AS{ty^ 



1/2 



. 1/2 

^ ) # log ^ 

T2 / oil ° T2 



\ 

satisfies the commuting equation^"^^ 

'I 



dL_ 

dtk 



■■ J 

(5.2.12) 

(5.2.13) 



(iii) The functions Tn{t) satisfy the following hilmear identity, forn > m+1, 
and all t,t' G C°° , where one integrates along a small circle about oo, 

T^{t - [z-^])T^+i{t' + [z-^])e^'^^'^-'>' z^-'^-Hz = 0. (5.2.14) 

(iv) The KP -hierarchy for = 0, 1, 2, . . . and for all n = 1,2 .. ., 



(_9_ld_l_d_ ^ 1 , 
^'^'^dh' 2 8^' 3dh' ■ • " 2dt,dtk+3 ) ° - 



of which the first equation reads: 



Remark: In order to connect with classical integrable theory, notice that, 
when r satisfies the equation above, the function 

q{ti,t2,...) := 2- 



dtl 



^^in terms of the Lie algebra splitting (j5.2.3|) . 

^•^ Remember the Hirota symbol (|5.2.4p and the Schur polynomial notation (|5.2.5p . 
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satisfies the classical Kadomtsev-Petviashvili equation (KP equation): 

d'^q d ( dq d^q dq\ r^oTa\ 

If q happens to be independent of t2, then q satisfies the Korteweg-de Vries 
equation 

dq d^q ^ „ dq /c: o 1v^ 



Proof: Identity (15.2.91) follows from the general fact that the product of two 
matrices can be expressed as a symmetric sum of determinant^, in particular 
the square of a Vandermonde can be expressed as a sum of determinants: 



e+k-2 

(fc) 



A\ui,...,Un) = det (^M^- 

Indeed, 

nhnit) = n\ det ninit) 



l<kl<n 



l<kl<n 



^ det ( / ^ Pt{za{k))dz„{k)] 
Yl / '^''^{^-tt''),^u.^ Ptiz^{k))dz^{k) 



^l{z)Ylpt{zk)dzk. 



k=i 



(i) At first, note 



Jr ^ z / 

/i»,i(t) + + ^/iij+2(t) + • • • , 

^Indeed, E<Tes„ ^^et (a,_^(j) ^j,^0))i<.,j<„ = ^^t (a^fc)i<,^fc<„ det (&^fc)i<,^fc<„ • 
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which by formula (15.1.81) of section 15.11 leads at once to the following rep- 
resentation for the monic orthogonal polynomials z) and their Stieltjes 
transforms, 

p.(uz) = ." """-If'" 

rn{t) 

M^l^^du = ,-nI^,±l(t±K!Jl. (5.2.18) 

Z-U Tn[t) 

(ii) The matrix L := SAS*"^ satisfies the standard Toda lattice. One com- 
putes 



— — = /ii+fcj implymg — — = A ^ 
Otk Oil 



Ik ULk 

Then, using the factorization (15.2.21) and the definition (15.1.91) of L = SAS~^ 
of section 15. H one computes 



L^ + —S-' + S 



Upon taking the ()_ and ()o parts of this equation {A_ means the lower- 
triangular part of the matrix A, including the diagonal and Aq the diagonal 
part) leads to 



Upon observing that for any symmetric matrix the following holds, 

ac\_/a ^\_r)f(^c\ f a c 
c h ) \ 2c h ) ^ \c h ) ~\c b 

it follows that the matrices L(t), S(t) and the vector p(t; z) = {pnit; z))n>o = 
S(t)x{z) satisfy the (commuting) differential equations and the eigenvalue 
problem 

dS 1 

— = --{L'),S, L{t)p{t- z) = zpit- z), (5.2.19) 
103 



and thus the tridiagonal matrix L satisfies the standard Toda lattice equa- 
tions 



dtk dtk 



^S-'SAS-'-SAS-'^S-' 



dt 



dtk 



with p{t; z) satisfying 

dp dS 



X{z) 



dtk dtk^''' 2' '"■"'^'■-^ 2' 

The two leading terms of pn(t; z) look as follows, upon using (15.1. 6p and 
(I5XT8|) : 



Pn{t;z) 



—^Pnit; z) 

Tn+l 

rn{t-[z-']) 



(5.2.20) 



'n+1 



z^'W-z 



^^dTrJdti 



+ ...]. 



(5.2.21) 



Thus, z"' admits the following representation in terms of the orthonormal 
polynomials Pi. 



p„ + ^^Z|^^«-i + 0(z"-2^ 



n+l 



dTn/dti 2\ 
Pn + Pn-1 + 0{Z ) 



(5.2.22) 



Then, using (15.2.211) in zpn and then using the representation (15.2.221) for 2;" 
and 2;"'"'"^, one checks that the diagonal entries 6„ and non-diagonal entries 
an of L are given by 

dTn/dti ■ 



hn = {zPn.Pn) 



Tn 



(z"+\p„)-(^",p„)- 



dTn+l/dti dTn/dti 



Tn 



— log 

Oti Tn 
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an = {zPn,Pn+l) = + • • • ,Pn+l) = 



and 



establishing (I5.2.12|) . 

(iii) The bilinear identity: The functions r„(t) satisfy the following identity, 
for n > m + 1, t,t' & C°°, where one integrates along a small circle about oo, 

Tn{t - [z-^])Trn+i{t' + [z-^])e^'^^'^-'>\''-'^-^dz = 0. (5.2.23) 

z=co 

Indeed, using the r-function representation for the monic orthogonal poly- 
nomials and their Stieltjes transform (15.2. ISp . one checks: 

rn{t)Tm{t ) A=oo 

^ / ^n '^njt- [Z "^]) ^-m T-m+i(f + [z^^]) ^y^oc dz 
Jz=oo ^n{t) Tmit') Z 

= i rf^eS?°(*-*^)^p„(t;z) /^^^^^eS?°*>>(«)rf« 

Jz=oo Jr Z — U 

= 2ni / e^^''^'~'''^'''pn{t]z)prn{t']z)e^^^'''''p{z)dz, using Lemma.^^ 
Jr 

= 2Txi I Pn{t] z)pm{t'; z)e^^ p{z)dz = 0, when m < n — 1, 
Jr 

by orthogonality, establishing (I5.2.23p . 

(iv) The KP hierarchy: Setting n = m + 1 in (15.2. 14p . shifting t t — y,t' ^ 
t + y, evaluating the residue, Taylor expanding in y^ (see (I5.2.6P ) and using 
the notation 

~ _ / 9 19 Id 
one computes the following residue about z = oo. 



= dze-^^^y^^\n{t-y-[z-'])rn{t + y+[z-']) 

= ^ 4,dze-^^ 2y,.'gEi — s^gEi y^a^rJt - u)rJt + 
2m J 



u 



u=0 
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OO 

e^i y^w^ J2 Sii-2y)si+iid)Tn o r„ 



OO 



CO / d d \\ 

for arbitrary y^, implying 

d / 92 ~ \ 

—r o r = and I ^^^^^ - 2Sk+i{d) 1 r„ o r„ = for = 1, 2, . . . . 

Taking into account the fact that trivially {d/dti)T o r = and that the 
equation above is trivial for k = 1 and k = 2, one is led to the KP hierarchy: 

In particular, for /c = 0, one computes 

1 1 

S4(t) := ^ + 2^2^!^ + tsti + 2^2^ + ^4, 

leading to the first equation in the hierarchy 



d \ ^ f dV 92 \ f d"^ V 



This ends the proof of Theorem I5.2[ ■ 

Remark: As mentioned earlier this method is very robust and can be general- 
ized to other integrals, besides (15.2.91) . upon using multiple orthogonal poly- 
nomials. Such integrals with appropriate multiple time-deformations lead 
to r-functions for multi-component KP hierachies; see Adler-van Moerbeke- 
Vanhaecke [121. 
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6 Virasoro constraints 



6.1 Virasoro constraints for /^-integrals 

Consider weights p{z)dz = e~^^^^dz with rational logarithmic derivative and 
E a disjoint union of intervals: 

-P- = V'{z) = § = and E = [\ [c2._i, c^,] C F C M, 

where F = [A, B] is an interval such that 

lim fiz)p{z)z^ = for all A; > 0. 

Consider an integral Init,c;P), generalizing fl5.2.9l) . where t := (ti,t2,...) 



and c = (ci,C2, ...,C2r); namely with a Vandermondc^ to the power 2/? > 
instead of a square, and omitting the n! appearing in (15.2.91) . 

» n 

/„(t,c;/?):= / |A„(;2)|2m(e^^°*»^^p(z,.)^^^fc) for n>0. (6.1.1) 

fe=i ^ ^ 

Then the following Theorem holds: 

Theorem 6.1 (Adler-van Moerbeke [S]) The multiple integrals LJt) '■= 
In{t,c;P) with Iq = 1, satisfy the following Virasoro constraints^ for all 
k > -1: 

E^^V(q)^ + E 41n(^,^) - h jKm.n(^>^)) m = 0, 

1 i>0 / 

(6.1.2) 

where ^'^^ ^i^il'^''"') ^'^^ combined differential and multiplication 

(linear) operators. For all n G Z, the operators 3i^^\{t,n) and Jfc^^(t,Ti) form 



25A„(z) = ^l<^<,<„(^^-^J) 



<i<.j<.n \ 

^^When i? equals the whole range F, then the the first term, containing the partials 
with respect to the Ci's, are absent in the formulae (6.1.3). 
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a Virasoro and a Heisenberg algebra respectively, interacting as follows 



'^k,ni '"£,11 

it(i) 



with "central charge" 



(2) 



12 



2/3 



1 - 6 iP'/' - p-'/'Y and c' = iQ-l). 



(6.1.3) 



(6.1.4) 



Remark 1: The operators J^^]^ = j[,^^(t, '^)'s are defined as follows: (the nor- 
mal ordering symbol ": :" means: always pull differentiation to the right, 
ignoring commutation rules) 

4!i = E ^ 4n jS : + (1 - /3) ((^ + 1) j£ - «S) , (6-1-5) 

i+j=k 

in terms of the — I^^\{t, n)'s. Componentwise, we have 

4!i(^, ^) = 4'^ + ^4°^ and ifl = = n^o, 

and 

jgl(t, n) = /3jf + {2n(5 + {k + l){l- (5)) J« + n{n(5 + l-[5) jf\ 
where 



7(0) _ 



5, 



kO 



J, 



(2) 



d 1 . ,^ 

i+j=k ■' —i+j=k ■' —i—j=k 



(6.1.6) 



The integer n appears explicitly in Jjf]j(i, n) to indicate the explicit n-dependence 
of the nth component, besides t. 
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Remark 2: In the case f3 = 1, the Virasoro generators (6.1.6) take on a 
particularly elegant form, namely for n > 0, 



lit) = 5^ : J«(t) : = Ji'\t) + 2nji'\t)+n'6, 

i+i=k 

4!iW = J^\t)+n5^ 



Ok 



JQk, 



wit 



i+j=k * ~i+j=k 



Ed 1 . . 



-i-j=k 



(6.1.7) 



One now establishes the following Lemma: 
Lemma 6.2 Setting 



n 

dinix) := |A„(x)|2^n (^^^*''^^-/'(^'^)^^'^ 
fc=l 



i/ie following variational formula holds: 



~;~dlji{xi I >■ + £ f {xi)Xi^ J 
de 



£ = 



oo 

5^ (a.4'i,n - rf/n. (6.1.8) 

1=0 



Proof: Upon setting 



£{x,t) ■.= X[e^^'^^lp{x,), 



(6.1.9) 



the following two relations hold: 



IE 



— Ok,0 



2 ^ dtidti 2 

ij>0 



9 



— + n5kfi S 
dtk 



+ ^ E 

l<a<l3<n l<a<n 



( ^ x^j^, allA;>0. (6.1.10) 

\l<a<n / 



^^The expression J^^"* = for fc = 
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So, the point now is to compute the e-derivative 

which consists of three contributions: 
Contribution 1: 



(6.1.11) 



£ = 



1 



d 



2/3|A(x)|2/5 9£ 



|A(x + £/(x)x'=+i)|"' 



e=0 



d 



1<Q<7<?1 



£ = 



l<a<7<n 
oo 

i=0 l<a<^<n 
/ 

oo 



£=0 



3y X '■y 



— 1) 



i+j=e+k 

i,j>0 
l<CK<'y<n 



l<a<n 



(9^ 



, 2 ^ dtidt. 2 

ij>0 



+ n 



k+£+l\ f d 



dtk^ 



+ n5k+e,o 



n{n — l) 



5k+i,o ]£, using f l6.1.10p . 



2 ^ dtidti 



i,j>0 



n 



k + i + l\ d n{n-l) 



dt 



(6.1.12) 



Contribution 2 : Using /(x) = 
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d 



£=0 



J2 if'Mxt' + {k + l[dxi 
1 1 

oo n n 

e=o a=l 1 

£-^J2i^ + k + l)ael-^ + n6k+i,o]£ll dx^. (6.1.13) 



Contribution 3 : again using /(x) = OjX* 



d 



oo n 



a=l 



i=l a=l 
oo n 



£=0 



0=1 i=l a=l 

^ 00 n n ^ 



v 



=0 a=l 



e>0 a=l 

i>l 



J 



~E^M 7)7 +^4+£+i,o)+E"^E 



9 



(6.1.14) 

As mentioned, for knowing (16.1.81) . we must add up the three contributions 
1, 2 and 3, resulting in: 

d 



de 



dln{xi Xi + ef{xi)x^^ ) 



(6.1.15) 



e=0 



E [(^-^w + (2^/5 + (^ + A; + 1)(1 - /3))4^, + n((n - l)[3 + 1)5, 
00 \ 



£=0 



Finally, the use of fl6.1.6p ends the proof of Lemma 16.21 
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Proof of Theorem 16'. J 1 The change of integration variable Xi ^ Xi + 
ef{xi)x^^^ in the integral (16.1. ip leaves the integral invariant, but it induces 
a change of limits of integration, given by the inverse of the map above; 
namely the q's in E = |J^[c2i-i, C2i], get mapped as follows 

qk.c,-£/(q)c(=+i + 0(£2). 

Therefore, setting 

r 

1 

we find, using Lemma 16.21 and the fundamental theorem of calculus. 



ending the proof of Theorem 16. 1[ 



6.2 Examples 

These examples are taken from [H El [7l [89]; for the Laguerre ensemble, see 
also HH] and for the Jacobi ensemble, see 1491. 



Example 1 (GUE) 

Here we pick 



p{z) = e-^(^) = e-^ , V = g/f = 2z, 
aQ = 1, 6o = 0, bi = 2, and all other Oj, bi = 0. 



Define the differential operators 



d 
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in terms of the end points of the set E = [jf^ [c2i_i, C2i] C M. From Theorem 
16.11 the integrals 

« n 

J„= / A„(z)2r[e-^'+^-i*'^^d^fc (6.2.2) 
satisfy the Virasoro constraints 

- Bkin = (-Jgi + 2tf]2,n) k = -1,0,1, .... (6.2.3) 

The first three constraints have the following form, upon setting F = log J„ 
(this will turn out to be more convenient in the applications), 



i>2 

PI \ ^ 



n 



i>l 



For later use, take linear combinations such that each expression contains 



-B-iF =: V^iF 



-IBoF =: V,F= 



2 



Example 2 (Laguerre ensemble) 

Here, the weight is 



yielding 




i>2 * V 




1 V- d \ ^ 
--> iU^\F 

i>i V 


~ y 


1 . d 

2 i>i ^^i+i 


2^ 2 



(6.2.4) 



/ z 
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Oo = 0, Oi = 1, 60 = ~CL^ bi = 1, and all other a,, bi = 0. 
Here define the (slightly different) differential operators 

Thus from Theorem 16.11 the integrals 

~ n 

4= / A„(^)2j]z,"e-^^+SS.*.4rf^, (6.2.6) 
"^s" k=i 

satisfy the Virasoro constraints, for k > —1, 

- Bkin = (-Jgl,n - «jSl,n + 4+2,n) (6.2.7) 

Written out, the first three have the form, again upon setting F = log/„. 




Replacing the operators Bi by linear combinations 

V2 = -Bo - {2n + a)B.i 

= -Bi- {2n + a)Bo- {2n + a)^B^i (6.2.8) 

yields expressions, each containing a pure derivative dF/ dti 
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V,F = —-^^u—-nin + a) 

OF ( \ 

Notice the non-linearity in this expression is due to the fact that one uses 
F — login rather than 7„. 

Example 3 (Jacobi ensemble) 

The weight is given by 

, . . V /I , x6 T./ 9 _ a-b+{a + b)z 
pab[z) e ^ [1- z) [1 + z) , V ^ fTr^2 

ao = 1, ai = 0, a2 = —1, 6o = a — 6, 6i = a + 6, and all other Oj, 6^ = 0. 
Here define 

The integrals 



/ An{zfl[{l-Zkni + Zkfe^^=^'*'Uzk (6.2.10) 
fe=i 

satisfy the Virasoro constraints {k > —1): 

- B,In = (4!),,„ - jgi + 6ojgi,„ + &ijg2,n) 4- (6-2.11) 
Introducing the following notation 

a — 2n + bi 
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the first four liaving tlie following form, upon setting F = log /, 



^ i>l ^^'^ i>2 * V 

^ ^ id , d \ - , d d . d"^ \ ^ 



9 9 



92 \ 9F(9F 
+2 F + 2 



dti dh dt2 dti+2' 



2-^\ F (f—Y-f—V 2——\ 



(6.2.12) 



7 Random matrices 



This whole section [7] is a very standard chapter of random matrix theory. 
Most of the results can be found e.g., in Mehta [69], Deift [32] and others. 

7.1 Haar measure on the space Tin of Hermitian ma- 
trices 

Consider the most naive measure (Haar measure) 

n 

dM:=Y[dMii Y[ d'^eMij d'^mMij (7.1.1) 

1 l<j<J<n 

on the space of Hermitian matrices 

Tin := {n X n matrices such that = M}. 
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The parameters Mij in (17.1.11) are precisely the free ones in M. This measure 
turns out to be invariant under conjugation by unitary matrices: (see Mehta 
|69| and more recently Deift [32] ) 



Proposition 7.1 For a fixed U G SU{n), the map 

Hn^nn-.Mi — >M' = U-^MU 

has the property 



dM = dM' I.e. 



det 



DM 



Proof: Setting M' = U'^MU, we have 

Tr = Tr M'^ 



and so 



Working out this identity, one finds 



1 i<i<i<" 

n 

1 l<i<j<n 



(7.1.2) 



Setting 

M := {Mii,...,Mnn,^e Mi2,...,3fJe M„_i,„, 3m Mi2,...,3m M„ 



identity fl7.1.2p can be written, in terms of the usual inner product ( , ) in 

n(n+l) 



(M, DM) = (M', DM') 
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for the x diagonal matrix with n I's and n{n — 1) 2's, 



/ 1 



D 



\ 



O 



o 



V 



2 / 



Let V be the matrix transforming the vectors M into M' 

M' = VM 

and so, fl7.1.2p reads 

{M',DM') = {VM,DVM) = {M,V^DVM), 
from which it follows that D = DV and so 



implying 



^ det D = detiV^DV) = (det Vf det D, 



detVl = 1. 



But for a linear transformation M' = VM the Jacobian of the map is V itself 
and so 

detV = 1, 



det 



dM 

ending the proof of Proposition 17.11 



Proposition 7.2 The diagonalization M = UzU^^ leads to "polar" or "spec- 
tral" coordinates M ^ {U, z), where z = diag{zi, . . . , Zn), Zi G M. In these 
new coordinates 

dM = A'^{z)dzi ...dzn dU. (7.1.3) 



Proof: Every matrix M G 7i„ can be diagonalized, 

M = UzU~\ 
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with U = e^ze~^ G SU{n) and 
^ = X] {akeicke - eek) + i huicki + e^)) G with at/, bk,e G 



l<k<e<n 



Then, using the definition of the measure dM, as in (17.1.11) . and using the 
fact that [A, z] is a Hermitian, one computes 



dM 



d(e^^e-^) = diz+[A,z] + 0{A')) 



A=0 



A=0 



Y\_dzi Y\ d ^G[A,z]jk d '^m[A,z]jk 

1 l<j<k<n 
n 

1 l<j<k<n ^<j<k<n 



A=0 



A=0 



Sinc(Q 

[A, z] = ^ (afc^efcf - e^, z] + i fefe^efei- + eft, 2;]) 



l<fc,£<n 



= ^ (ofcK^^ - ^fe)(efe£ + e^fc) + 2 - ^fe)(efe£ - e^fc)), 

l<fc/<n 

with 

3fte[v4, z\k/ = Ukiize - Zk), Qm[A, z]ke = bki{z(, - Zk), 

estabhshing flTX^ near M = z. By Lemma [7H dM = d{U'~^MU') for any 
unitary matrix U\ implying the result (17.1.31) holds everywhere, establishing 
Lemma 17.21 ■ 

Remark: The set of Hermitian matrices is the tangent space to a symmetric 
space G/ K = SL{n, C)/ SU{n). The argument in Theorem 17.21 can be gener- 
alized to many other symmetric spaces, as worked out in van Moerbeke 
page 324-329. 



^Cki denotes the matrix having a 1 at the entry (fc,£) and everywhere else. 
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7.2 Random Hermit ian ensemble 

Consider the probability distribution on the space Hn of Hermitian matrices, 
in terms of Haar measure dM, given by 

P{M e dM) = J-e-^^WdM. 

Let Zi < Z2 < . . . < Zn he the real eigenvalues of M. Then 

P{zi e dZi, . . . ,Zn e dZn) = P{zi, . . . , Zn)dZi . . . dZn 

^A'{z)fle-^^^^^dz, 



Z, 



i=l 

with 



' Zl<...<Z„ 

Lemma 7.3 



Zn^ P{Z1, . . .,Zn)YldZi. 

J Zl<...<Zr,. 1 



» n 

P(M e Tin, spectrum (M) C ^) = i . (7.2.1) 



p It, 



Proof: Indeed 

P(M e Tin, spectrum MCE) 



1 /■ " 

= — / Al{z^,...,z^)llxAz^)e-''^'^^dz, 

J zl<...<Zn j^^l 

If " 

" 7r65""'Mi)<-<Mn) i=l 
1 /• " 
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^ 1 

n 

1 



showing Lemma 17751 ■ 

Let pq{z) , pi{z) , p2{z) , ... be orthonormal polynomials with regard to the 
weight p{z) defined on R, as discussed in section 15.11 



Pi{z)pj{z)p{z)dz = 6ij, 
and Pn{z) be the monic orthogonal polynomials, 

Pi{z)pj{z)p{z)dz = hiSij. 



Then, 

Pn{z) = ^{Z" + ...) = -^Pn{z). (7.2.2) 

Proposition 7.4 Setting 

» n 
Zn= I Al{z)Wp{Zk)dZk, 

we have the identity 

Z-^Al{z)\\p{zk) = -^det{Kn{zk, ^^))i<M<n, (7.2.3) 
1 

where the symmetric kernel Kn is given by 

n~l 

Kn{y,z) = \/ p{y)p{z) ^Pj{y)Pj{z) (Christoffel-Darboux) 



j=0 



= ./^y/M^^"^^^^"'^^'^ '^"'^^'^^"^'^ (7.2.4) 
V hn-1 y- z 

The kernel Kn{y , z) has the following "reproducing" property 

/ Kn{y,u)Kn{u,z)du = Kn{y,z) and / Kn{z,z)dz = n. 
Jm. Jk 
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Proof: Notice that the Vandermonde An{z) can also be expressed as 
dct{pi_i{zj))i<ij^n by row operations, where Pi{z) can be chosen to be any 
monic polynomial of degree i, and in particular the monic orthogonal poly- 
nomial of degree i. Thus, one computes for the normalization Z^, 

» n 
Zn ^ A\z)l[p{Zi)dZi 

n 

det(pj_i(^j))i<jj<„ det{pk-i{ze))i<k,e<nYlp{zi)dzi 

i=l 

n „ 

= (~l)''"^'''n / PAk)-ii^k)P7r'{k)-liZk)p{Zk)dZk 

n,7r'eSn k=l "^^ 

n „ 

= "^Yi Pl{k)-\{zk)p{zk)dzk , using the orthogonality of the p^'s, 

7re5„ fc=i 

n—l p n—l 

= n! JJ / pl{z)p{z)dz ^ nlYlhk- 



Then using the expression obtained for Z„ and (det74)^ = det (AA^) in the 
third equality, one further computes 

n 

Z-'Al{z)l[piz,) 

1 

1 " 

= . T-rn-l , (Pi- 1 ) ) 1<» det {pk-l{Zi)) l<k/<n TT P(^fc) 



^et feM,/;^ det fciMyMi^ 



l<i,j<n \ V ^ / i<fe,£<n 



— det{Kn{Zk, Zi))i<k/<n- 
77/. 



l<A;,^<n 
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One finally needs the classical Christoffel-Darboux identity: setting pj (y) = 
for J < 0, one checks 

n—l n—l 

iy - z) ^Pj{y)pj{z) = {ajj_ipj^i(y) + ajjPj{y) + ajj+iPj+i{y)) Pj{z) 
j=o 

n-l 
j=0 

= On-l.n iPniy)Pn~l{z) - Pn-l{y)Pn{z)) , 

and one uses (17.2.21) . The reproducing property follows immediately from 
the Christoffel-Darboux representation of the kernel in terms of orthogonal 
polynomials. This proves Proposition 17. 4[ ■ 



7.3 Reproducing kernels 

Lemma 7.5 Let K{x,y) be a symmetric kernel satisfying the reproducing 
property 

/ K{x,y)K{y,z)dy = K{x,z). 

JR 

Then 

j det {K{zi,Zj))^^^ .^^ dzn 

= i^j K{z,z)dz -n+l^det{K{zi,Zj))i<ij<n-i, (7.3.1) 
where dy stands for any measure on M. 

Proof: The proof of (17.3.11) is due to J. Verbaarschot [90], which proceeds in 

two steps: 

Step 1 : Let 



Mn 



Mn-1 m 



m7 7 
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he a n X n Hermitian matrix, with M„_i an — Ixn — 1 Hermitian matrix. 
Then m is a (n — 1) x 1 column and 7 G M. Then 

det M„ = 7 det M„_i — m^M„_im. 

Indeed given the column u G C"^^, one checks 



/ W Mn-i m \ I I u 
1 I [ Th^ 7 j I 1 



Mn-l Mn-lU + m 

VL^ Mn-l + ffl^ uJ Mn-lU + ffl^ U + VL^ TTl + 7 



Mn-l 

7 - m^M, 



^1 ^ ^ , upon setting m = — Af^i^^m. 



The latter assumes that M„_i is invertible. Furthermore, using Af^i^^ 
(detM„_i)-iAf„_i, 



det Af„ 



(7 - m ' M-l^m) det Af, 



n— 1 



(7 - m^Af„_im(det Af„_i)-^) det Af„_i 
7 det A^n-i — frJ Mn-im. 



Step 2 : Proof of identity (17.3.ip . Define the matrix 



and 



y K{zk-i,Zk) ) 
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Mk is a symmetric matrix, since is a symmetric kernel. One finds, upon 
integration over M, 

AQt{Mn)dZn 

« « n— 1 

= det(M„_i) / K{z,z)dz- I dZn K{Zn, Zi){Mn-l)ijK{Zj, Zn) 

» n— 1 „ 

= det(M„_i) / K{z,z)dz -'^{Mn^l)ij i K{Zn, Zi)K{Zj, Zn)dZn 

« n—l 

= det(M„_i) / ^(M„_i)ij(M„_i),-, 

» n—l 

= det(M„_i) / K{z,z)dz-det{Mn-i)y2(^n-i)ij{Mn-i)ji 

= det(M„_i) (^J K{z, z)dz - {n - 1)^ , 
establishing Lemma [7.51 ■ 
Lemma 7.6 // 

{i) K{x,y)K{y,z)dy = K{x,z) 



(ii) / K{z, z)dz = n, 
Jr 

then 

. det{K{zi, Zj))i<ij<ndzm+i .. .dzn = {n- m)\ det(i^'(2:i, -2j))i<ij<m 

(7.3.2) 

Proof: The proof proceeds by induction on m. On the one hand, assuming 
(17.3.21) to be true, integrate with regard to Zm and use identity (I7.3.ip : 



/ . . . / det{K{Zi, Zj))i<ij<ndz.mdZr, 

J Jr"-™+i 
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(n-m)! / dzmdet{K{zi,Zj))i<ij<m 

{n-m)\{ K{z, z)dz - m + 1) Aei{K{zi, 2i))i<jj<-m-i 



= {n- m)\{n - m + 1) dei{K{zi, 2;j))i<j,j<m-i 
= (n - m + 1)! det{K{zi, %))i<ij<m-i- 

On the other hand for m + 1 = n, the statement follows at once from Lemma 
17. 51 This ends the proof of Lemma 17.61 □ 

7.4 Correlations and Fredholm determinants 

For this section, see M. L. Mehta [69], P. Deift [32] , Tracy- Widom [83| and 
others. Returning now to the probability distribution on the space Tin of 
Hermitian matrices (setting p{z) := e"^*-^-*) 

P(M e dM) = ^e-^'^^^^dM, 

remember from Lemma 17.31 and Proposition 17.41 



P(M e Hn, spectrum M C E) = / Pniz)dzi ...dzn 

with 

n 

Al{z)Y[p{z,)dzi 

Pn{z) = - = Z;,'Al{z)X{p{z,) 

^l{z)\{p{z,)dz, 
1 

= —<lei{Kn{zk,zt))i<k,t<n, (7.4.1) 
n\ ~ ~ 

with the kernel Kn{y,z) defined in Proposition 17. 4[ 

n-l 

Kn{x,y) = ^p{x)p{y)^pj{x)pj{y). (7.4.2) 

j=0 
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Let E be the expectation associated with the probabihty P above. Then one 
has the following "classical" Proposition for any subset E dM. (for which a 
precise statement and proof was given by P. Deift [32]): 

Proposition 7.7 The 1- and 2-point correlations have the following mean- 
ing: 

/ Kn{z, z)dz = E,{^ of eigenvalues in E) 
Je 

/ det{Kn{zi, Zj))i<ij<2dzidz2 = of pairs of eigenvalues in E), 

Jexe 

and thus ^ 

Kn{z,z) = — E (t^ of eigenvalues in dz) (7.4.3) 

Proof: Using (17.3.21) for m = 1 and (I7.2.3p . one computes 

Kn{z,z)dz 



XEizi)Kn{zi,Zi)dzi 

/ dziXsizi) ■■■ det{Kn{zi, Zj))i<ij<ndz2 ...dz„ 

JR J JR"-l 



(n-1 



[ dziXsizi) [ ■■■ [ ■^Al{z)f\p{zk)dz2...dz„ 

n 

Xe(2:i)A2(z) Jjp(zfc)c/Zfc 

1 

If " 

— / #{z such that Zi e E}Al{z)T\p{zk)dzk 

jRn ^ 



n 
n 



E of eigenvalues in E) 



29 



If xi,X2 G E, then it counts for 2 in the second formula. 
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A similar argument holds for the second identity of Proposition I7.7[ ■ 

Consider disjoint intervals Ei, . . . , and integers 1 < ni, . . . , < n 
and set n^+i := n — nj. Then for the nxn Hermitian ensemble with P„ 
as in (I7.4.1I) . one has: 

P (exactly rii eigenvalues G -Ej, 1 < i < m) 



n 

Til, •••) ^m? 



\ „ ni 711+712 
^ / -^K" i=l i=ni+l 



711 + . ..+?!„ 



i=ni + ...+7i™,_i+l i=J]™»ifc+l 

(7.4.4) 

This follows from the symmetry of Pn{x) under the permutation group; the 
multinomial coefficient takes into account the number of ways the event oc- 
curs. 

Lemma 7.8 The following identity holds 

^ n / m \ 

k=l V 7=1 ' / 

m 

I + Kn{x, y)^XiXE^ (y) 



det 



7 = 1 



Proof: Upon setting 



K{x,y) := Kn{x,y)^XiXESy)^ 



and upon using the fact that K{x, y) has rank n in view of its special form 
(17.4.21) . the Fredholm determinant can be computed as follows: 

det{I + k{x,y)) 



SZf! /^^det[i^(xi,Xj 



l<i,j<i 



dxi...dxe 



KiyXi, Xj^ 



dxi...dxi 
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dxi...dxi 



k=l 

£=0 ■ l<si,...,Sf<m "^'^^ r=l 

» n 

det[Kn{xi,Xj)]i<ij<ndxi...dxn, using Lemma [731 



=0 l<si, ...,S£<m 
1 



?!(n-£)! 



£=0 l<si,...,Sf <ni l<ii<...<i£<n 

P„(x)(ixi...(ix„ 



„ n / m \ 

/ n + X -^^^s (^'^^ Pn{x)dXi...dXn, 



establishing Lemma [7.81 ■ 

Proposition 7.9 T/ie Fredholm determinant is a generating function for the 
probabilities: 

P{exactly rii eigenvalues G -Ej, 1 < i < m) 



ra, ! V dX 



I + ^\iKn{x,y)xE^{y) 



.(7.4.5) 



a// Ai=-i 



In particular 

P{no eigenvalues G Ei, 1 < i < m) = det / — y)^^™^;. (y) • (7.4.6) 

Proof: The first equahty below follows from Lemma 17.81 Concerning the sec- 
ond equality below, in order to carry out the differentiation YYi ;n ( ^ ) ' 
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one chooses (keeping in mind the usual product rule of differentiation) a first 
group of Hi factors, a (distinct) second group of n2 factors, ... , a mth group 
of rim factors and finally the last group of n — rii — . . . — rim remaining factors 
among the product 11^=1 + YlT=i ^iXs i^k)^ ■ Then one differentiates the 
first m groups, leaving untouched the last group, where one sets Aj = — 1. 
This explains the second equality below. Let C"^ „^ be the set of distinct 
committees of size ni,n2, ■■■,nm, n^+i '■= n — rii — ...—rim formed with people 



1,2, ...,n: 



nA \dXi 



I + ^XiKn{x,y)xE^ (y) 



all A,=-i 



i=l vc^i/ JR"k=l \ i=l J 

„ ni ni+...+nm 



all Ai=-i 



i=ni+...+nm~i 



Y[ I ^ ~ 5Z ] Pn{x) dXi...dXn 

i=m+...+nm+l \ i=l 

( n \ ^ 



rii, ...,nm,n- y^Ui 

V 1 J 



rti 



ni + ...+n„ 



i=l 



i=ni+...+nm-i 



JJ ( 1 - y^^XeS^i) ) Pn{x)dXi...dXn 

i=ni+...+n,„+l \ £=1 / 

P(exactly rii eigenvalues G Ei, 1 < i < m), 



as follows from 07.4.41) . thus establishing identity fl7.4.5p . whereas (17.4.61) 
follows from setting ni = . . . = Um = 0, completing the proof of Proposition 
El ■ 
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8 The distribution of Hermitian matrix en- 
sembles 



8.1 Classical Hermitian matrix ensembles 
8.1.1 The Gaussian Hermitian matrix ensemble (GUE) 

Let Tin be the Hermitian ensembles 

'Hn = {n X 77, matrices M satisfying = M}. 

The real and imaginary parts of the entries Mij of the n x Hermitian matrix 
(M = M""") are all independent and Gaussian; the variables Ma, 1 < i < n, 
3fJe Mij and Mij, 1 < i < j < n, which parametrize the full matrix have 
the following distribution (set ?i„ 3 H = [uij), with real Ua and off-diagonal 
elements Uij := Vij + iwij) 

1 _ 2 

P(Mii e duii) = —=e ^"duii 1 <i <n 

'71 



2 2 

P(3fJe Mjk e dvjk) = —=e~'^''^''dvjk I < j < k < n 

'71 



2 2 

P(Q=m Mjk e dwjk) = —j=e~'^'^:>>'dwjk I < j < k < n. 

'TT 



Hence, using Haar measure (17. 1.11) . 
P(M G dH) 

n 

= Y[P{M,, e duii) Yl P{^e Mjk e dvjk)P{Qm Mjk e dw,k) 

1 i<j<k<n 

n n 

1 l<j<k<n 1 l<j<k<n 

n 

1 l<j<fc<n 

= Cne-^''"dH = CnAl{z)l[e-'"dz„ HeHn, (8.1.1) 
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using the representation of Haar measure in terms of spectral variables Zi 
(see Proposition I7.2p and where 




This constant can be computed by representing the integral of (18.1.11) over 
the full range as a determinant of a moment matrix (as in (15.2.91) ) of 
Gaussian integrals. 



8.1.2 Estimating covariances of complex Gaussian populations 
and the Laguerre Hermitian ensemble 

Consider the complex Gaussian population x = (xi, . . . , Xp)^ , with mean and 
covariance matrix given by 

= E(x) = (/ii, . . . ,/ip)^, S = {E{xi - fii){xj - fij))^<ij<p 

and density (for the complex inner-product ( , )), 

1 ^^l[g-ft,S-^x-fi)) 



-e 



(27r)p/2(detS)V2 

Let Ai > A2 > . . . > Ap > be the eigenvalues of E. Taking n samples of 
{xi, . . . , Xp)~^ , consider the normalized p x n sample matrix: 

/ 2:11 - a^ii) Xi2-^{^^Xu) ... 

X — 

\ Xpi — - {'^i Xpi^ Xp2 — - {'^i Xpij . . . Xpn — - Xpi) / 

and the p x p sample covariance matrix, 
1 

1' 

which is an unbiased estimator of S. It is a classical result that when S = /, 
the eigenvalues 2:1, . . . , 2;^ > of S" have the Wishart distribution, a special 
case of the Laguerre Hermitian ensemble (see Hotelling [S3] and also Muir- 
head M) 



S := — — -xx'^ , with eigenvalues zi, Zp > 0, 



P„,p(5 e dM) = CnpAl{z)Y[e-'^z'l-''-'dzidU = e-^'^\det My-^-^dM. 
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8.1.3 Estimating the canonical correlations between two Gaussian 
populations and the Jacobi Hermitian ensemble 



In testing the statistical independence of two complex Gaussian populations, 
one needs to know the distribution of canonical correlation coefficients. I 
present here the case of real Gaussian populations, not knowing whether the 
complex case has been worked out, although it should proceed in the same 
way. To set up the problem, consider p + q normally distributed random vari- 
ables {Xi, Xp)~^ and (Yi, ...,Yg)~^ {p < q) with mean zero and covariance 
matrix 



The method proposed by Hotelling [SJj is to find linear transformations U = 
LiX and V = L2Y of X and Y having the property that the correlation 
between the first components Ui and Vi of the vectors U and V is maximal 
subject to the condition that Var Ui = Var Vi = 1; moreover, one requires 
the second components U2 and V2 to have maximal correlation subjected to 



etc. . . 

Then there exist orthogonal matrices Op E 0{p), Oq E 0{q) such that 



p 




Ip 




(i) Vart/2 = Var^a = 1 

(ii) U2 and V2 are uncorrected with both Ui and Vi 



J:^|^'^^2^22^' = O^PO, 



where P has the following form: 
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Pk 



o 



Pk+l 



i y 



o 



o 



\ 



Pp 



p, k = rank S 



12, 



1 > Pi > P2 > • • • > Pfc > 0, pk+1 = ... = Pp = (canonical correlation coefficients), 

Pi are solutions (> 0) of det(S]j^^Si2S^2^E72 — p^I) = 0. 
Then the covariance matrix of the vectors 



U = LiX := OpS^i^/^X and V = L2Y := O^T^-^'^Y 



-1/2, 



has the canonical form (det Scan = Hill ^ Pi)) 
cov I , , 1 = S 

with 



spectrum 
and inverse 



1, . . . , 1, 1 - pi, 1 + pi, . . . , 1 - Pp, 1 + Pp 

q-p 



^-1 



1 



-P 



From here on, we may take S = Scan- The n (ri > p + g) independent 
samples (xn, . . . , Xip, yn, . . . , yig)'^, . . ., (a;„i, . . . , x„p, ?/ni, • • • , J/ng)"^, arising 

from observing ^ y- ^ l^ad to a matrix ^ ^ of size (p + g, ri), having the 

normal distribution [70] (p. 79 and p. 539) 



(27r)-«(p+'?)/2(detS)-"/2exp -Ixr 



^11 ^12 

Sj^2 ^22 



(^27]-)-"(P+9)/2(det S)""/2e"^ (xT(E-l)nx+yT(E-i)22y+2yT(s-l)T x) 
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The conditional distribution oi p x n matrix x given the q x n matrix y is 
also normal: 

(det27rf])-"/2e-5^""'(^-^^)(^-^^)^ (8.1.2) 

with 

Q = Sii - Ei2Sj2^E2i = diag(l - pj, . . . ,1 - p^) 

P — Ei2E22^. 

Then the maximum likelihood estimates of the pi satisfy the determinantal 
equation 

det{S^^^Si2S^2^Sj2 - r^I) = 0, (8.1.3) 

corresponding to 



12 "^22 



^ _ , 5*11 5*12 \ ( xy 



SJ2 S22 J ' V y^^ yy^ 



where Sij are the associated submatrices of the sample covariancc matrix S. 

Remark: The r j can also be viewed as = cos 6i , where the 61, ...,6p are the 
critical angles between two planes in M": 

(i) a p-dimensional plane = span {{xn, Xni), {xip, ...,Xnp)} 

(ii) a g-dimensional plane = span {{yn, ...,y„i)'^, {yig, ...,y„q)}. 

Since the (g,n)-matrix y has rank(?/) = q, there exists a matrix Hn G 0{n) 
such that yHn = {yi O); therefore acting on x with Hn leads to 

q n—q q n—q 

yi/ = (yi I O) I xH^ = {u\v)\p . (8.1.4) 

With this in mind, 

'S'l2'S'22^'S'2'2 — r^'S'ii 

= xy^{yy^)~^yx^ —r'^xx^ 

= xH{yHy {yH{yHy)-^yH{xHy - r\xH)(xHy 

T\\-l /T\ /T 





){'• 


I, 


' 







uu 





I ^) ( ^ ] { It ] -r'{u I 



M ' \ 2f I \ u 

'u \ v) \ 
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and so the equation (18.1.31) for the rj can be rewritten 

det{uu^ -r^{uu^ + vv'^) = 0. (8.1.5) 

Then setting the forms (I8.1.4P of x and y in the conditional distribution 
(18. 1.21) of X given y, one computes the following, setting H := Hn, 

TtQ-\x - Py){x - Pyy 

= Tin'\xH ~ PyH){xH - PyHy 

= Tin-' {{u I v) - P{yi I O)) {{u I v) - P{yi \ O))^ 

= Trn-\u- Pyi)iu- Py^y + TTn-'vv'^; = diag(l - p^, . . . , i _ 

this establishes the independence of the normal distributions u and v, given 
the matrix y, with 

u = NiPy,,n), v = N{0,n). P = diag(pi,...,pp). 

Hence uu~^ and vv~^ are conditionally independent and both Wishart dis- 
tributed; to be precise: 

• The p X p matrices vv'^ are Wishart distributed, given y, with n — q 
degrees of freedom and covariance f2; 

• The p X p matrices uu'^ are non-centrally Wishart distributed, given y, 
with q degrees of freedom, with covariance Q and with non-centrality 
matrix 

^Py,yJP^n-\ 

• The marginal distribution of the q x q matrices yy~^ are Wishart dis- 
tributed, with n degrees of freedom and covariance Ig, because the 
marginal distribution of y is normal with covariance Ig. 

To summarize, given the matrix y, the sample canonical correlation coef- 
ficients rf > . . . > are the roots of 

(r^ > . . . > Tp) = roots of det{xy^ {yy^)^^yx^ — r^xx^) = 
= roots of det(MM^ — r'^{uu^ + vv^)) = 
= roots of det (uu^ (uu^ + vv~^)~^ — r'^I) = 0. 
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Then one shows that, knowing uu^ and vv'^ are Wishart and conditi- 
ionally independent, the conditional distribution of rf > . . . > r^, given the 
matrix y is given by 

1 

^ (g/2).C.(/,) |A|! ^^(^ 

wherj^ 

i?^ = diag(ri, . . . ,rp), c 



2 _ j;„„^„2 „2n ^ rp(?7-/2) 

' '^'^'^ r,(g/2)r,((n-g)/2)r,(p/2)' 

and where the C\ are proportional to Jack polynomials corresponding to the 
partition A; for details see Muirhead [70] and Adler-van Moerbeke [11]. By 
taking the expectation with regard to y or, what is the same, by integrating 
over the matrix yy~^ , which is Wishart distributed (see section 18. 1.21) . one 
obtains: 

Theorem 8.1 Let Xi, . . . , Xp, Yi, . . . ,Yq (p < q) be normally distributed 
random variables with zero means and covariance matrix Tj = ( 

V ^21 ^22 

If Pi, ... , pp are the roots of det(S]j^^Si2S22^^72 ~ P^^) = 0; ^^^'^ maxi- 
mum likelihood estimates , . . . , from a sample of size n (n > p + q) are 
given by the roots of 

(let{xy^ {yy^Y^yx^ — r^xx'^) = 0. 

Then, assuming pi = . . . = p^ = 0, the joint density of the Zi = rf is given 
by the following density: 

T,^'l\^,p,,/\p{z)\{z^~^-'^'\l - z,Y-^-^-'^/'dz,. (8.1.6) 



^"^Using the standard notation, defined for a partition A, 

{a)\ := ]^(a + (1 - , with (a;)„ := x{x + 1) . . . (a; + n - 1), xq ^ 1. 
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Remark: Taking complex Gaussian populations should introduce in the for- 
mula above Ap{z) instead of Ap{z) and should remove the 1/2's in the ex- 
ponent. 

8.2 The probability for the classical Hermitian random 
ensembles and PDE's generalizing Painleve 

8.2.1 The Gaussian ensemble (GUE) 

This section deals with the Gaussian Hermitian matrix ensemble, discussed 
in previous section. Given the disjoint union of intervals 

r 

E:=[j [c2i-i,C2i] CM, 
1 

define the algebra of differential operators 

B.-Y.^% (8.2.1) 

i=l 

The PDE ( I8.2.2P appearing below was obtained by Adler-Shiota-van Moer- 
beke [H El l^, whereas the ODE fl8.2.3l) was first obtained by Tracy- Widom 
|84j . The method used here is different from the one of Tracy- Widom, who 
use the method proposed by Jimbo-Miwa-Mori-Sato [5l]- John Harnad then 
shows in [51] the relationship between the PDE's obtained by Tracy- Widom 
and by Adler-van Moerbeke. 

Theorem 8.2 The log of the probability 

„ n 

/ Al{z)lle-^?dz. 

P„ := P„(a// ZieE) = 1 

/ Al{z)\{e-^hz, 

satisfies the PDE 

{Bti + 8nB\ + 12Bl + 2ABo-l6B^iBi)\ogFn + Q{B\\ogFr,y = 0. (8.2.2) 
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In particular 

f{x) := — logP„(max2i < x) 
ax i 

satisfies the 3rd order ODE 

f" + 6 /'2 + 4{2n - x^)f' + 4xf = 0, (8.2.3) 
which can be transformed into the Painleve IV equation. 

Proof: In Theorem 15. 2[ it was shown that integral (here we indicate the t- 
and i?-dependence) 

1 r " 

r^{t;E) = - Al{z)r\e~'?+^^'''^'dz, (8.2.4) 
satisfies the KP equation, regardless of E, 



(8.2.5) 

and in Theorem 16. 11 Tn{t; E) and r„(t; M) were shown to satisfy the same Vira- 
soro constraints, and, in particular for the Gaussian case, the three equations 
(16.2.41) . with the boundary term missing in the case r„(t; R). 

Let Tj denote the pure t-differentiation appearing on the right hand side 
of fl6.2.4p . with "principal symbol" ^ 



dU+2 ' 

d 1 V- . d 



1 • 



dti 2^ dti^i 

^ i>2 * ^ 



d 1^ .. d 
Ti := -K- 



d 1 \ - d n \ ^ d 

r. yiti- yiti^ • 8.2.6 

dts 2^ 'dti+i 2^ 'dti.i ^ ' 
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Recall the differential operators Vi in terms of the boundary operators fl8.2.ip . 
appearing in the Virasoro constraints fl6.2.4p . 



= -\b, 

V, = (8.2.7) 

With this notation, the Virasoro constraints f l6.2.4p can be summarized as 
{F:=\ogTn{t-E)) 

V.,F = T^,F-^, VoF = ToF~—, V,F = T,F-^. 
Expressing the action of Tj on ti, one finds 

T_iti = l-t2 T\t^=T_,{l-t2) = ^t3 

T^ti = -nt2 T\ti = T_^T^_^ti = r_i(^t3) = -3t4, 

one computes consecutive powers of "Dj and their products, and one notices 
that Vi involves differentiation with regard to the boundary terms only, im- 
plying in particular that Vi and Tj commute. In view of the form of the KP 
equation, containing only certain partials, and in view of the fact that the 
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''principal symbor of Tj equals d/dti+2, it suffices to compute 
rati 

= T_iF 

2 

r, / nil \ r, n , 



nil 
nil 



ViF = TiF 



V_{DiF = V^iTiF = TiV_iF = Ti [ T_iF 

-.2 



nHi 



T^^F 1, 



T\F + —U. 



n 



VoF = ToF - 



VlF = VoToF = ToVoF = To TqF - 



n 



T^F. 



(8.2.8) 



Since one is actually interested in the integral fl8.2.4p along the locus C : = 
{all = 0}, and since readily from (18^.61) one has Ti\c = d/dti^2, one 
deduces from the equations above 



V\F 
V\F 



VlF 



V_iViF 
VoF 



T\F 

t\f 



n d'^F 



T^F 



c 2 dtl 

- A 
c dtl 

- A 
c dt2 



T\F 
ToF 



n 

d 



d^F 



c V dt2 



c dt\ 
F 



TiT_iF 
ToF 



= ^T_iF 
c dts 



c \dt3dt1 2dt2 



_ dF 

£~y ~ ^2 



n 



^^Notice one also needs VqF, because dFjdt-i appears in the expressions V^F and 
V^iViF below. 
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By solving the five expressions above hnearly in terms of the left hand side, 
one deduces 



dtl 
dF 

d^F 



V\F 



VnF 



n 

2 



n 



V_iViF 



d^F 



dt\ 
d'^F 



dtl 



V\F 



VlF 



+VoF 



So, substituting into the KP equation and expressing the Vi in terms of the 
Bi as in (18 .2 .71) . one finds 



n 



c 



+ 6 Vi^F 



n 



V\F + 6nV\F + 3VlF - SVqF - AV^iViF + 6(V\Fy 



^{B\ + 8nB\ + 12Bl + 24Bo - 16B^iB,)F + l{B\Ff 
io o 



which establishes f l8.2.2p for (remember notation f l8.2.4p ) 

F = logr„(0,E) = log P„ (all Zi G + log r„(0, M). 

Since the Bk are derivations with regard to the boundary points of the set E 
and since logr„(0,M) is independent of those points, the equation f l8.2.2p is 
also valid for logP„; it is an equation of order 4. 

When is a semi-infinite interval (— oo,x), then one has Bk = x'^^^d/dx 
and then, of course, PDE (18.2.21) turns into an ODE (I8.2.3p . of an order 
one less, involving /(x) := ^ log P„ (maxj < x). For the connection with 
Painleve IV, see section 18. 3[ thus ending the proof of Theorem |HI 
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8.2.2 The Laguerre ensemble 

Given E C M"*" and the boundary operators 

^'^^=E^^'a^' for fc = -1,0,1,..., 

the following statement holds: (see [H El [7] for the PDE obtained below; the 
ODE was first obtained by Tracy- Widom [8^) 

Theorem 8.3 The log of the probability 

n. n 

P„ := P„(a// z^eE) = ^- 

/ Al{z)l[zte-^'dz, 

satisfies the PDE 



f B\ - 2i33 ^ + (1 - o?)B\ - ABiB^i + ?>Bl 
\ +2(2n + a)i3oi3_i - 2Bi - {2n + a)i3o 




+ Q{B\ logP„)2 - A{B\ logP„)(^-i fogPn) = 0. 

(8.2.9) 

In particular, /(x) := log P„(maXj 2;j < x) satisfies 
x'^f'" + xf' + 6xf^-4ff'-{{a-xf-Anx)f-{2n + a-x)f = 0, (8.2.10) 
which can be transformed into the Painleve V equation. 

8.2.3 The Jacobi ensemble 

The Jacobi weightn is given by (1 — z)"-{l + z)^. For E C [— 1,+1], the 
boundary differential operators Bk are now defined by 

Introduce the following parameters 

r = A{a^ + l?), s = 2{a^-b'^), q = 2{2n + a + bf. 
r = a^ + b'^, s = a^-b'^, q={2n + a + bf. 
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Theorem 8.4 (Maine- Semengue and Adler-van Moerbeke ^) The fol- 
lowing probability 

„ n 

/ A^izfUil-z^ril + z^fdz, 
P„ := P„(a// G = ^4 (8.2.11) 

/ A^{zf\[{l-z,Y{l + z;)'dz, 
■^[-1.1]" k=i 



satisfies the PDF: 

Bt^ + (g - 2r + 2)B^, + q{?>Bl - 2Bo + 2B2) + AB^Bl^ 
-2{2q - l)BiB-i + (2i3_ilogP„ - s){Bi - + 2BoB-i] 



logPn 



+2B\ log P„ {2Bo log P, + 3B\ log P„) = (8.2. 12) 

In particular, f{x) := (1 — x^)^ logP„(maXj Aj < x) for < x < 1 satisfies: 

{x^ - Iff" + 2{x^ - 1) {xf" - 3/'') + (8x/ - q{x^ - 1) - 2sx - 2r) f 

-f{2f-qx-s) = 0, (8.2.13) 

which is a version o/Painleve VI. 

Proof of Theorems \8.3\ and \8.4\ ' It goes along the same lines as Theorem 18.21 
for GUE, namely using the Virasoro constraints fl6.2.8p and (16.2.91) . together 
with the KP equation (KTB . This then leads to the PDE's (KTM and 
08.2. 12p . The ODE's (18.2.101) and fl8.2.13|) are found by simple computation. 
For connections with the Painleve equations see section 18.31 ■ 



8.3 Chazy and Painleve equations 

Each of these three equations (l8X3|) . (18.2.101) . 08.2.131) is of the Chazy form 

r + + p'^ - '-^ff + ^/^ + - ^/ + ^ = 0, (8.3.1) 

with P,Q,R having the form: 

Gauss P{x) = 1 4Q(x) = -Ax^ + 8n R = 

Laguerre P{x) = x 4Q(x) = — (x — af + 4nx R = 

Jacobi P(x) = l-x2 4Q(x) = -(^(x^ - 1) + 2sx + 2r) R = 

144 



The differential equation fl8.3.ip belongs to the general Chazy class 

= F{z, f, /', /"), where F is rational in /, /', /" and locally analytic in z, 

subjected to the requirement that the general solution be free of movable 
branch points; the latter is a branch point whose location depends on the inte- 
gration constants. In his classification Chazy found thirteen cases, the first of 
which is given by equation (18.3.11) . with arbitrary polynomials P{z), Q{z), R{z) 
of degree 3,2,1 respectively. Cosgrove and Scoufis |30l |29], show that this 
third order equation has a first integral, which is second order in / and 
quadratic in /", 

f" +^ {(Pf" + Qf + R)f' - (P'f" + Q'f + R')f 

+ liP"f' + Q")f - lP'"f + = 0; (8.3.2) 

c is the integration constant. Equations of the general form 

= Gix, /, /') 

are invariant under the map 

aiZ + 02 , „ 05/ + + 07 
X ^— s> and J I— s> . 

a^z + 04 a^z + 04 

Using this map, the polynomial P{z) can be normalized to 

P{z) = z{z — 1), z, or 1. 

Equation (18.3.21) is a master Painleve equation, containing the 6 Painleve 
equations, replacing f{z) by some new variable g{z). e.g, 

• g"^ = -Ag'^ - 2g'{zg' -g) + Ai (Painleve II) 

• ^"2 ^ _4^/3 ^ 4(^^/ _ gY ^ ^^^/ ^ (Painleve IV) 

. {zg"f = {zg' - g) {-Ag'^ + A,{zg' - g) + A^) + A,g' + A4 

(Painleve V) 
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. {z{z - l)g"y = izg' - g) (Ag'^ - Ag\zg' -g) + A^) + A^g'^ + A,g' + A, 



Now, each of these Painleve II, IV, V, VI equations can be transformed into 
the standard Painleve equations, which are all differential equations of the 
form 



whose general solution has no movable critical points. Painleve showed that 
this requirement leads to 50 types of equations, six of which cannot be re- 
duced to known equations. 

9 Large Hermitian matrix ensembles 

9.1 Equilibrium measure for GUE and Wigner's semi- 



Remember according to (17.4.31) . the average density of eigenvalues is given by 
Kn{z, z)dz. Pastur and Marcenko [76] have proposed a method to compute 
the average density of eigenvalues (equilibrium distribution), when n gets 
very large. For a rigorous and very general approach, see Johansson [55] . 
who also studies the fluctuations of the linear statistics of the eigenvalues 
about the equilibrium distribution. 

Consider the case of a random Hermitian ensemble with probability de- 
fined by 



for a diagonal matrix A = ( ). Consider then the spectral function 

of A, namely da{\) := - 5(A — Oj). The Pastur-Marcenko method tells 
us that the Stieltjes transform of the equilibrium measure of dv{\), when 
n ^ oo, namely 



(Painleve VI) 



/" = F{z^ /, /'), rational in /, /', analytic in z. 



circle 





satisfies the integral equation 




(9.1.1) 
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The density of the equihbrium distribution is then given by 

— — = -Qm f{z). 

dZ IT 

When A = 0, the integral equation flQ.l.ip becomes 

1 



-z - v'^f{z) 



with solution f(z) = ^{—z± \/ z^ — 4t>2), and thus one finds the classical 
semi-circle law, 

1 



dz TT 



53m f{z) 



for l^;! > 2v, 



concentrated on the interval [—2v,2v]. 

As an exercise, consider now the case, where v = 1 and where the diagonal 
matrix A has two distinct eigenvalues, namely 

A = diag(a, . . . , a, /3, . , /? ) . 

pn {l-p)n 

See e.g., Adler-van Moerbeke [I3]. The integral equation (19.1.11) becomes 

/- ^ ^"^ , ^ = 0' forO<p<l, 

(j- z- f a- z- f 

which, upon clearing, leads to a cubic equation foi g := f + z, 

- {z + a + P) + {z{a + P) + aP + 1) g - a^z - (1 - p)a - pP = 0, 

having, as one checks, a quartic discriminant Di (z) in z. Since the roots of a 
cubic polynomial involve, in particular, the square root of the discriminant, 
the solution g{z) of the cubic will have a non-zero imaginary part, if and only 
if Di{z) < 0. Thus one finds the following equilibrium density, 

^jy^^-^ ^ J —Qm g{z) for z such that Di{z) < 



dz TT 



'^m f{z) 



for z such that Di{z) > 
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Therefore the boundary of the support of the equihbrium measure will be 
given by the real roots of Di{z) = 0. Depending on the values of the pa- 
rameters a, (3 and p, there will be four real roots or two real roots, with a 
critical situation where there are three real roots, i.e., when two of the four 
real ones collide. The critical situation occurs exactly when the discriminant 
D2{a,P,p) (with regard to z) of Di{z) vanishes, namely when 



D2{a, P,p) = 4p(l - p)j{j^ - + 37(9p2 _ Qp + 1) _ i)S 



= 0. 

7=(a-/3)2 



This polynomial has a positive root 7, which can be given explicitly, the 
others being imaginary, and one checks that, when one has the relationship 

a — (5 — — =, upon using the parametrization p 



^, ^j^v^ii ^^^^^^ ^^^^ i^^^^^^^^^^^^^K,^^^^ ^ ■? I 1 ' 

V^g^ - g + 1 r+1 

two of the four roots of D-i{z) coUide. This is to say, this is the precise 
point at which the support of the equilibrium measure goes from two to one 
interval. This then occurs exactly at value 

on the real line. 

9.2 Soft edge scaling limit for GUE and the Tracy- 
Widom distribution 

Consider the probability measure on eigenvalues Zi of the n x n Gaussian 
Hermitian ensemble (GUE) 



P„(all Zi&E)^ 



n 

^2 



' Al{z)l[e-^fdz, 



n 

^2 



Al{z)l[e-^Uz, 



Given the disjoint union E := lj^[a;2t-i, X2i\ C M, define the gradient and the 
Euler operator with respect to the boundary points of the set E: 

2f Q d 
^- = Ea^ S^^Y."^,—. (9.2.1) 
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Remember the definition of tlie Fredliolm determinant of a continuous kernel 
K{x,y), the continuous analogue of the discrete kernel fl2.U.3p . 



(—1)" f 

det(/ - K{x, y) XE{y)) = 1 + J] / ^-^i ' ' ' det(ir(zi, z^)) 

n=l 

We now State: [Ell [H [5] 
Theorem 9.1 The gradient 

f{Xu...,X2r) := V.l0gP(E'=), 

with 



F{E^) := lim P„ (all V2n^ (zi - V2^] G E' 

n—>oo \ \ / 

satisfies the 3rd order non-linear PDF: 



Vi - A{S, -2^)f + 6(V./)' = 0. (9.2.2) 
In particular, for E = {x, oo), 

J^{x) := lim P„ (V2n^ f^^nax - v^) < x) 
= det(/ - A X(:,,oo)) 

= exp ^— J {a — x)g^{a)da^ , (9.2.3) 
is the Tracy- Widom distribution, wzi/j^ 



, A{x)A'{y)-A'{x)A{y) f°°^^ 

A{x,y):=—-^ — — ^ ^ = / A{u + x)A{u + y)du 

x-y Jo 

and g{a) the Hastings- McLeod (unique) solution of 

g" = ag + 2g^ 

_2^l (Painleve II). 

~ ^W7? fora / oo. 



■^^ Remember the Airy function: 



A{x) 

TT 



1 f°° fu^ \ 

— I COS [- xu\ du, satisfying the ODE A"{x) = xA{x). 

^ Jo V 3 / 



149 



Proof: Step 1: Applying Proposition 17.9^ it follows that 

P„(all Zi G E') = det (J - Kn{y, z)xAz)) 



where the kernel Kn{y, z) is given by Proposition 17.41 



1(2^2+^2) Pn{y)Pn~l{z) - Pn-l{y)Pn{z) 

2) y-z 



(9.2.4) 



(9.2.5) 



The p„'s are ortho normal polynomials with respect to e ^ and thus propor- 
tional to the classical Hermite polynomials: 

1 „ . . 1 



Pniy) 



2"/2V^7rV4 



Hn{y) 



--y" + 



(9.2.6) 



with 



Step 2: The Plancherel-Rotach asymptotic formula (see Szego [H2]) says that 



2n/2+l/4^^!7rl/4 



x=^/2n+T+ 



A{t) + Oin~^/'), 



V2 nl/6 



uniformly for t G compact i^' C C and thus, in view of (19.2.6^ . 



e 2 Pn{X} 



2^/%-^2 (^A{t) + 0(n"i)) . (9.2.7) 



Since the Hermite kernel (19.2.51) also involves p„_i(x), one needs an estimate 
like (19.2.71) . with the same scaling but for p„ replaced by Pn-i- So, one needs 
the following: 

t 



X 



V2n + 1 + 



(2n-l)(l + - -) + 



t 



V2n-1 ( 1 + ^ ^ , + O' ^ 



2n-r v^(n- 1)1/6(1 + _i_)i/6 

t 



V2n - 1 + 



2n- 1 
v^(n - 1)1/6 



O 



+ 



7/6 / ■ 



v^(n - 1)1/6 



1+ 0( 



n 
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Hence, from fl9.2.7l) it follows that 



e 2p„_i(xj 



e 2 p^^^i^x) 



V2 nl/6 



v^{n-l)l/6 



(9.2.8) 

From the definition of the Fredholm determinant, one needs to find the limit 
of Kniy] z)dz. Therefore, in view of (19.2.51) and using the estimates (19.2.71) 
and dnXSD, 



lim Kniy; z)dz 



-lim iiy\-'^^y'^^'^ 



Pn{y){Pn{z) -Pn-l{z)) - Pn{z){Pn{y) - Pn-liv)) 



{y - z)V2 nV6 



ds 



lim f-V^'(2V4n-ii)2n-V3 

n^oo V 2 / 



A{t) {A{s + 72-1/3) - A{s)) - A{s) {A{t + n-i/3) - A{t)) 



n' 



-V3(t 



ds 



A{t)A\s) - A{s)A'{t) 



ds = A{t, s)ds 



t - s 

and thus for E := \jil^[x2i-i,X2i] and E := [Jl^[c2i-i,C2i], related by 

E 



E = V2^ + 



(9.2.9) 
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one has shown (upon taking the hmit term by term in the sum defining the 
Fredholm determinant) 



hm P„(all Zi e 



hm P„(all Zi e V2n + 



E' 



hm P„(aU G (2^ + 1)^/^ + 



E' 



hm det(/ - 

n— >oo 

det(/ - AXe)- 



£;=(2ii+l)l/2 + 



Step 3: From Theorem 18.21 P„(an Zi G E'^) satisfies, with regard to the 
boundary points Cj of E, the PDF fl8.2.2p : thus setting that scahng into this 
PDF yields (remember Vc and are as in (19.2.11) and Bk = '^i'^^'^r. 



= {B\ + SnB'i^ + 12Bl + 2ABo - IQB^iBi) \ogFn 

+ 6(i3^1ogP,,)' 



Ci=V2n+ 



1, 



+ o(n2/3)_ 



Note that in this computation, the terms of order n^^^ cancel, because the 
leading term in 



12Bl - I6B-1B1 



+ 



cancels versus the leading term in SnB'^^ = 16n^/^V^ + . . .; thus only the 
terms of order n^^^ remain. Since in step 2 it was shown that the limit 
exists, the term in brackets vanishes, showing that logP(£''^) satisfies the 
PDF ( 19:2:21) . 

Step 4- In particular, upon picking E = (x, 00), the PDF (19.2.20 for 
fix) = ^log ^"(0;) = ^log hm P„ (y2n^ (^z^ax - V2nj < 



dx 

becomes an ODF: 



r-4xr+2/+6r = 0. 
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Multiplying this equation with /" and integrating from — oo to x lead to 
the differential equation (the nature of the solution shows the integration 
constant must vanish) 

r' + 4/'(r-x/' + /) = 0. (9.2.10) 

Then, setting 



j2 ^^2 ^4 (9.2.11) 



/' = -g' 

f = g'^ - xg^ - g" 

and, since then /" = ~2gg' , an elementary computation shows that (19.2.111) 
is an obvious solution to equation fl9.2.10p . For (19.2. lip to be valid, the 
derivative of the right hand side of the second expression in (19.2.110 must 
equal the derivative of the right hand side of the first expression in (19. 2. lip , 
i.e. we must have: 

= (/)' -f = {g" - xg' - g'y + g' = 2g'{g" - 2g' - xg\ 

and so g" = 2g^ + xg. Hence 

— logJ^ix) = f' = -g^ 

Integrating once yields (assuming that g^ decays fast enough at oo, which 
will be apparent later) 

— \ogJ='{x) = J g^{u)du; 
integrating once more and further integrating by parts yield 

/»00 rOO /"OO ^ /"OO 

logJF(a;) = / da g'^{u)du = / da{^—a) I g'^{u)d' 

J X J a J X J a 



U 



CO 



oo 



oo 



= a g {u)du + ag {a)da 

J a ^ J X 

POO 

= / {a-x)g^{a)da, (9.2.12) 

J X 

confirming (19.2.30 . For x — > oo, one checks that, on the one hand, from the 
definition of the Fredholm determinant of A, the two leading terms are given 
by 

J^(x) = det (/ - AX^^ = [ '^)^^ + • • ■ 

J X 
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and, on the other hand, from fl9.2.12p . the two leading terms are 

/>oo 

J-'{x) = 1 — {a — x)g'^{a)da + . . . 

J X 

Therefore, comparing the two expressions above, 

{a — x)g'^{a)da = / dzA{z,z) + ... 



and upon taking two derivatives in x, 



-g'^ix) = ^A(x,x) + ... 



— / A{u + xfdu + ... 
d 

— J A{ufdu + ... 



= -AixY + ..., 

showing that asymptotically g{x) ~ A{x) for a; — >• oo. It is classically known 
that asymptotically 

_2 3/2 OO 

A[x) = — ;=— TT-ll + > , +...), as X ^ OO. 

The solution g{x) of g" = 2g^ + xg behaving like the Airy function at oo is 
unique (Hastings-McLeod solution). It behaves like 

g{x) = Ai(x) + O ( ^^^^ I for X — > oo 



X 



The Tracy- Widom distribution JF of mean and standard deviation (see Tracy- 
Widom [HZ]) 

E{J^) = -1.77109 and a{J^) = 0.9018 
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has a density decaying for x oo as (since J-'{x) tends to 1 for x oo) 

poo roo T 

/ / -L 4 3 /2 

J-''{x) = J^{x) I g'^{a)da ~ / A'^{a)da ~ - — e~3^ for a; — >• oo. 

Jx Jx OTTX 

(9.2.13) 

The last estimate is obtained by integration by parts: 

4^3/2 °° ^ 1 /-OO °° 

o /-, X — ^ ^7 \ r / c -J 



Then, foUowing conjectures by Dyson [35] and Widom [93], Deift-Its-Krasovsky- 
Zhou [33J and Baik-Buckingham-DiFranco [18] give a representation of J-'{x) 
as an integral from — oo to a; and thus this provides an estimate for a; — > — oo, 

H^) - 2./V<-')£Jpe.p{£(fi(,)-i,= + i-).,} 



where 

f oo 



POO 

R{y)= / g{sfds = g\yf-yg{yf-g{y)\ 

J V 
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